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The vertical intensities of mesons penetrating 5.25 and 30 cm of lead have been measured 
to an altitude of 15,000 feet and of those penetrating 20 cm of lead to an altitude of 32,000 
feet (275 millibars pressure) by coincidence counter telescopes sent up in an aeroplane from 
Bangalore, magnetic latitude 3.3°N. A comparison of our results with those of Schein, Jesse, 
and Wollan indicates that the latitude effect between 3.3°N and 52°N of the vertical intensity 
of mesons shows no marked increase even to altitudes corresponding to pressures of 275 milli- 
bars. This is in striking contrast with the total intensity, which shows a very pronounced 


increase of latitude effect to these heights. 


INTRODUCTION 


HE variation with altitude of the total 
vertical intensity of cosmic rays has been 
measured by several authors'~’ at different lati- 
tudes using coincidence counter telescopes. The 
results of such experiments have definitely estab- 
lished the existence of a large latitude effect at 
high altitudes and are in satisfactory agreement 
with the deductions of the cascade theory. For 
brevity we shall call any curve giving the vertical 
intensity of some component or group of cosmic 
rays as a function of altitude the altitude-in- 
tensity curve of that component or group. Valu- 
able information on the processes leading to the 
creation, annihilation, and scattering of mesons 
1 Pfotzer, Zeits. f. Physik 102, 23 (1936). 
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and their dependence on energy could also be 
obtained through measurements of the altitude- 
intensity curves of mesons of different energies. 
There are very few such measurements. Dymond® 
measured the altitude-intensity curve of mesons 
penetrating 10 cm of lead at Edinburgh, magnetic 
latitude 59°N, but the work was interrupted by 
the war and only a preliminary note has been 
published. Schein, Jesse, and Wollan,®!° and 
Schein, Wollan, and Grétzinger™ have carried 
out similar measurements at Chicago, magnetic 
latitude 52.5°N, using lead absorbers of different 
thicknesses, and the results of these authors are 
the only ones available yet. Since cosmic-ray 
mesons are all produced in the atmosphere, a 
comparison of the altitude-intensity curves for 
mesons of different energies at different latitudes 
would give valuable information about the de- 

’ Dymond, Nature 144, 782 (1939). 

®M. Schein, W. P. Jesse, and E. O. Wollan, Phys. Rev. 
57, 847 (1940). 
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pendence of their creation on the energy of the 
primary rays. It is therefore important to know 
the altitude-intensity curves of mesons at differ- 
ent latitudes, and particularly near the geomag- 
netic equator. No such measurements have been 
made previously near the geomagnetic equator. 

As is well known, India is particularly suited 
for such measurements. First, the geomagnetic 
equator passes through south India. Secondly, 
due to the fact that the magnetic field of the 
earth is not axially symmetric with regard to its 
geographic center, India lies in the region of the 
earth’s surface where the cosmic-ray intensity is 
a minimum, according to the world survey by 
Millikan and Neher." This means that only the 
most penetrating cosmic rays can reach the sur- 
face of the earth in south India. Even theoreti- 
cally (cf. Johnson") the critical energy required 
by a cosmic ray to reach the earth at the mag- 
netic equator in India is over 15 X 10° ev whereas 
on the equator west of India this critical energy 
is lower, and reaches the minimum value of 
13 10° ev in South America. Theoretically the 
critical energy should reach its maximum value 
in some region in the Pacific Ocean at about 
140°E longitude, that is north of Australia. 
However, the earth’s magnetic field is not that 
of a pure dipole and it is possible that owing to 
local variations in the field south India lies in 


- the region where the critical energy is a maxi- 


mum, as the sea-level survey of Millikan and 
Neher seems to indicate. In any case, south 
India lies sufficiently near, if not in, the region 
of maximum critical energy for a particular 
importance to attach to measurements made 
here of the altitude-intensity curves of mesons 
near the geomagnetic equator. 

Vertical intensity- measurements if properly 
carried out are capable of attaining an accuracy 
sufficient for a quantitative comparison with 
theory. It is desirable for a proper comparison 
with theory that in experiments designed to 
measure the vertical intensity the geometry of 
the apparatus should be such that the lengths of 
paths traversed in the atmosphere and the 
absorber by the different rays measured by it do 
not differ by more than about 10 percent. Al- 


2R. A. Millikan and H. V. Neher, Phys. Rev. 47, 207 


(1935). 
13 T, H. Johnson, Rev. Mod. Phys. 10, 193 (1938). 
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though such a condition tends to reduce the 
counts per minute, this should be remedied by 
increasing the number of counters in a tray, and, 
where possible, by making measurements for a 
relatively longer period, rather than by reducing 
the distance between the extreme counters and 
thus spoiling the geometry, and consequently 
the accuracy of the measurements. Pfotzer'’s! 
arrangement for the measurement of the total 
vertical intensity satisfies this condition but the 
apparatus used by other workers does not always 
do so. For example, in the experiments of Neher 
and Pickering* for the measurement of the total 
vertical intensity, the maximum allowed angle 


was 47° from the vertical entailing a path length ~ 


47 percent greater than that of a vertical particle. 
Their geometry was further complicated by 
putting two counters in parallel at a horizontal 
distance roughly equal to half of the vertical 
separation of the counters. : 

It has been shown by Greisen and Nereson“ 
that the effect of side showers on counter tele- 
scope measurements can be considerable. Since 
the vertical meson intensity increases much more 
slowly with altitude than that of the soft com- 
ponent, it is always necessary in measurements of 
meson intensity to reduce the effect of side 
showers to a minimum. According to Greisen and 
Nereson the effect of side showers for a triple 
coincidence counter telescope is about 14 percent 
at an altitude 4300 meters and it would be greater 
at higher altitudes. A quadruple coincidence 
counter telescope is therefore desirable in meas- 
urements of vertical meson intensity which at 
great heights is only a small fraction of the 

tal intensity. 

In order to obtain the maximum information 
possible from the various experiments it is 
necessary that the results of different authors 
should be easily comparable and reproducible. 
For example, quantitative information about the 
latitude effect of the different components: can 
only be obtained if this is possible. To this end 
we suggest that in executing and communicating 
all future experiments for measuring the varia- 
tion with altitude of the vertical intensity of 
cosmic rays of different penetrating power certain 
criteria be adopted which we indicate below. 


4K, Greisen and N. Nereson, Phys. Rev. 62, 316 (1942). 
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SCALE IN CM 
Fic. 1. Arrangement of counters. 


The cosmic-ray intensity at any point in the 
atmosphere is mainly a function of the amount of 
matter traversed by the rays in reaching that 
point. Hence comes the importance of measuring 
the vertical rays only and keeping the variation 
of path length within about 10 percent by 
suitably designing the geometry of the apparatus, 
as we have already mentioned above. Now the 
pressure is a direct measure of the mass of air 
traversed by a vertical ray, and it is therefore 
appropriate that the intensity should be given 
as a function of pressure rather than in meters 
of water equivalent, as is often done. It is also 
what is always measured directly in the experi- 
ment. The conversion of pressure into meters of 
water equivalent involves certain physical as- 
sumptions about the relative absorbing powers 
of air and water and the small uncertainty it 
introduces is undesirable, especially in view of 
the fact that the accuracy of cosmic-ray measure- 
ments will continually increase in the future. A 
measurement of pressure in centimeters of mer- 
cury is however arbitrary and inconvenient for 
theoretical purposes. We therefore propose that in 
the future the intensity be given as a function of 
pressure measured in millibars. Since 1000 milli- 
bars, which is equal to the pressure of the 
standard atmosphere, corresponds to between 10 
and 10.3 meters of water equivalent, the pressure 
in millibars divided by 100 gives immediately to 


within 3 percent the thickness of the intervening 
atmosphere in meters of water equivalent. This 
procedure would therefore have the advantage of 
communicating what is directly measured in any 
experiment in a form convenient both for a 
comparison with theory and with the past experi- 
mental results. 

It is not the absolute intensities, but the 
relative values of the intensities which are of 
significance in comparing the work of different 
authors. It is therefore essential that the reading 
at sea level both with and without any absorber 
that may be used should always be measured 
with the apparatus concerned and communicated 
in the paper. Since the latitude effect at sea 
level and the absorption curve in lead at sea level 
at different latitudes is known, this would enable 
a fairly quantitative comparison of the curves of 
different authors to be made. We would also 
suggest that the final readings should always be 
given in a table, even if a graph is included, and 
the method of averaging the readings and the 
statistical accuracy be indicated. 

We have planned a number of experiments on 
the basis of the foregoing discussion and carried 
out a preliminary measurement which is being 
reported here. In this preliminary experiment, 
we have measured the variation with height of 
the intensity of cosmic rays penetrating 20- and 
30-cm thick blocks of lead absorbers at Banga- 
lore, magnetic latitude 3.3°N, by sending up 
quadruple coincidence counter telescopes in an 
aeroplane. The intensity of mesons penetrating 
5.25 cm of lead was also measured by Bhabha’s 
method!*'* with a triple coincidence counter 
telescope. Two flights were made, one on De- 
cember 26, 1944 to an altitude of 15,000 feet and 
one on December 28, 1944 to an altitude of 
32,000 feet. Although further experiments are 
still in progress, we think that the preliminary 
results which have already been obtained are of 
sufficient importance to merit being communi- 
cated now. 


APPARATUS AND MEASUREMENTS 


The counters used in the experiment were 
filled by a slightly modified form of the procedure 


18 Bhabha, Proc. Ind. Acad. Sci. Al9, 23 (1944). 
16 Chandrashekhar 
(1944). 
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given by Neher" and had a plateau of nearly 300 
volts. They were 15 cm long and 1}’’-1}” in 
diameter. The distance between the extreme 
counters was in every case 35 cm. Therefore no 
single particle making an angle of more than 22° 
to the vertical could traverse the telescope. Thus 
a particle recorded at this maximum angle of 
22° would only travel a thickness of the atmos- 
phere and absorber 8 percent greater than a 
particle arriving vertically. 

Three distinct sets were used for the measure- 
ments. One was designed on the basis of Bhabha’s 
method with a 5.25-cm thick lead absorber which 
allowed the measurement of the intensity of 
mesons with energies above 0.93 X 108 ev. It had 
three counters in a vertical plane working in 
coincidence with two additional counters in 
parallel on each side of the middle counter and 
put in anti-coincidence with the previous three. 
A lead absorber 1.25 cm thick was inserted below 
the top counter and a 4.0-cm lead. block was put 
between the middle and bottom counters. A side 
view of the arrangement is shown to scale in 
Fig. 1. It ‘has been shown" that this splitting of 
the lead absorber into two parts together with 
_ the anti-coincidence arrangement in fact allows 

the soft component to be cut out by only 5.25 cm 
of lead and that the two anti-counters to the 
sides of the middle counter are adequate for 
reducing the effect of side showers. 


OOO 
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SCALE IN CM 
Fic. 2. Counter telescope. 
17 op Stony an others; Procedures in Experimental Physics 


(Prentice-Hall, Inc., New York, 1938). 
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The second set had a quadruple coincidence 
counter telescope with 20 cm of lead as shown 
in Figs. 2 and 3. In the side view the horizontal 
lines indicate the active length of each counter, 
The third set was similar to the second one but 


-had, in addition to the 20 cm of lead, another 


10 cm of lead on the top. Each counter telescope, 
its amplifier, and its plate and bias batteries 
were all mounted in a box which was closed on 
all sides and waxed except for a small hole on 
the top made for the purpose of equalizing the 
pressure inside and outside the box. Inside each 
box was an electric heater operated from the 
aeroplane supply. The high voltage to the 
counters was supplied from a separate battery 
box which also had the necessary provision for 
heating. 

The pressure was measured by a standard 
altimeter. The outside temperature was also 
measured throughout the flight and was required 
in order to calculate the real altitude from the 
indicated altitude. The measurements were made 
during ascent and descent in both flights. In 
addition, the airplane flew for half an hour at 
each “indicated” altitude of 5000, 10,000, 15,000, 
20,000, 25,000, and 30,000 feet corresponding to 
pressures of 832, 671, 530, 430, 332, and 265 
millibars, respectively. 

The cosmic-ray counts of each of the three 
sets were recorded on three telephone call 
counters, and these together with a clock and 
the altimeter were photographed at the beginning 
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Fic. 3. Counter telescope. — 
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TaBLeE I. Readings taken on airplane flights. 


with om of lend endings with of tend Readingn with om of tend 
altitude Pressure in in counts in counts in 
in feet cm millibars counts min per hr. counts min. per hr counts min. counts/hr. 
2900 68.8 906 42 60 42.0+ 4.3 36 60 36.0+4.0 24 60 24.0+3.3 
5000 63.2 832 22 30 44.0+ 6.2 20 30 40.0+6.0 15 30 30.0+5.2 
7500 56.4 742 40 38 63.14 6.3 40 51 47.0+5.0 16 27 35.6+5.9 
10000 51.0 671 32 30 64.04 7.6 20 30 40.0+6.0 17 30 34.0+5.5 
12500 45.5 599 60 40 90.0+ 7.8 60 59 61.0+5.2 21 28 45.0+6.5 
15000 40.3 530 58 31 112.3410.2 36 31 70.0+6.7 22 31 42.6+6.1 
20000 32.7 430 81 47 103.4+7.6 
25000 25.2 332 146 58 148.2+8.2 
30000 20.9 275 274 71 231.5+8.6 


and end of each of these intervals, as also at 
intermediate points on the way up and down. 
The results on the way up and down were con- 
sistent. No heating was required up to 15,000 
feet and the heaters were not put in operation in 
the first flight. During the second flight two 
heaters failed to operate, and measurements up 
to 32,000 feet were obtained only of the intensity 
of mesons penetrating 20 cm of lead. The readings 
were consistent on the way up and down, and 


Pressure in millibors 


Fic. 4. Curve A—Vertical intensity of mesons pene- 
trating 30 cm of lead. Curve B—Vertical intensity of 
mesons penetrating 20 cm of lead. Curve C—Vertical in- 
tensity of mesons penetrating 5.25 cm of lead. (Note: Since 
the counters were not of identical diameter in all the three 
cases, quantitative comparison of the three curves together 
is not to be made.) 


in good agreement with those made during the 
first flight. The combined readings of the two 
flights are given in Table | and plotted in Fig. 4. 


INTENSITY OF MESONS OF DIFFERENT ENERGIES 


Certain definite conclusions can be drawn from 
our results given in Table I and Fig. 4. 

1. The intensity of cosmic rays capable of 
penetrating 30 cm of lead increases very slowly 
with altitude. The intensity at 15,000 feet 
(indicated) (530 millibars) is about twice the 
intensity at 2900 feet (905 millibars). 

2. The intensity of mesons capable of pene- 
trating 20 cm of lead also increases slowly with 
altitude, but at a more rapid rate than in 
case (1). 

3. The intensity of mesons capable of pene- 
trating 5.25 cm of lead, i.e., having an energy 
greater than 0.93 X 108 ev, increases perceptibly 
more rapidly with altitude. The intensity at 
15,000 feet is nearly 2} times the intensity at 
2900 feet. 

Schein, Jesse, and Wollan'® report that there 
is no difference in the intensity of rays pene- 
trating any absorber thickness between 4 and 
18 cm of lead and this appears to contradict the 
result of Schein, Wollan, and Grotzinger." We 
are unable to understand the reason for their 
having found no absorption of rays between 4 
and 18 cm of lead contrary to the theoretical 
expectation and our own experimental results 
and those of Bostick,!* Rossi and Greisen,!® and 
Hall** who find a great increase in the percentage 
of slow mesons with altitude. 


18 W. H. Bostick, Phys. Rev. 61, 557 (1942). 
19 B. Rossi and K. Greisen, Phys. Rev. 61, 121 (1942). 
2” D. B. Hall, Phys. Rev. 66, 321 (1944). 
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Pressure in millibors 


Fic. 5. Curve A: Vertical meson intensity at 3.3°N 
(Bhabha, Aiya, Hoteko and Saxena). Curve B: Vertical 
meson eee at 52.5°N (Schein, Jesse and Wollan— 
1941). Curve C: Vertical total intensity at 3.3°N (Neher 
and Pickering). Curve D: Vertical total intensity at 49°N 
(Pfotzer). : 


It is worth mentioning that our readings at 
10,000 feet for all the three apparatuses lie 
markedly below the three smooth curves drawn 
through the readings at the other altitudes, but 
there is not sufficient evidence to decide whether 
this result is purely a coincidence due to statis- 
tical fluctuations or is real. 


LATITUDE EFFECT 


In Fig. 5 we have plotted our results giving the 
intensity of mesons penetrating 20 cm of lead 
at 3.3°N as curve A and, for comparison, the 
latest results of Schein, Jesse, and Wollan!® for 
the intensity of mesons at 52.5°N as curve B. 
The two curves have been fitted at sea level to 
allow for the known latitude and longitude 
effect’? of 12 percent. In the same figure, we 
have plotted as curve C the variation of the 


total vertical intensity with altitude at 3.3°N as - 


given by Neher and Pickering*® for a triple 
coincidence counter telescope, this curve being 
fitted to our curve so as to show a ratio of vertical] 
meson intensity to total vertical intensity of 80 
percent as observed at ground level at Bangalore, 
Curve D gives the variation of the total vertical 
intensity with altitude’as measured by Pfotzer! 
at a magnetic latitude of 49°N. The four curves 
together show at a glance the striking fact that 
whereas the latitude effect between 3.3°N and 49°N 
of the total intensity shows a pronounced increase 
with altitude, the penetrating component shows 
practically no such increase of latitude effect even 
to heights corresponding to pressure of 275 milli- 
bars. The theory of Hamilton, Heitler, and Peng”! 
according to which the penetrating component 
should show only a slightly greater latitude 
effect than at sea level up to heights correspond- 
ing to a pressure of 100 millibars is therefore at 
least in qualitative agreement with our measure- 
ments. The difference in the geometry of the 
counter telescopes used by the different authors 
and the statistical accuracy of the results do 
not yet permit a quantitative comparison. 
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On the Scattering of Slow Mesons 
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Photographs of mesons of total energy lying between 1.55 10® and 2.55 X 108 ev have been 
obtained by a special arrangement of counters controlling the expansion of the cloud chamber. 
The scattering of these particles, both multiple and single, has been studied in two thicknesses 
of lead (2 cm and 4 cm). It is found that though the +/# law of Williams holds good for the 
average angle of scattering, its absolute value is only about 50 percent of the theoretical value 


expected from William’s formula. The Gaussian distribution of the number of particles with 


angle of scattering is found to be approximately true. A cross section for the non-Coulombian 
nuclear scattering has been calculated. The value of this is found to be 1.84 10~-* per nucleon 
for mesons of mean total energy 2X 10* ev. This agrees satisfactorily with that calculated by 
Bhabha and Weinberg and Ma for transversely polarized mesons of this energy, but does not 
agree with that found experimentally by Shutt and Code. Our value is about twenty-five 
times the value obtained by Shutt. The too low value of Shutt and others may be attributed 
to the presence of high energy mesons (for which the nuclear scattering cross section is very 
small) in large numbers so as to mask the effect of low energy mesons. One case of a high energy 


proton being singly scattered has also been obtained. 


INTRODUCTION 


XPERIMENTS on the scattering of meso- 
trons by plates of copper, lead, and gold 
have been done by Wilson! and also by Code,? the 
latter using a 3.5-cm tungsten plate as the 
scattering material. In these experiments the 
energy of the mesons photographed varies from 
2X10* to about 210° ev. Both these experi- 
ments confirm the Gaussian distribution of the 
number of particles as a function of the product 
of the energy E and the scattering angle @, and 
they also confirm within experimental errors the 
theoretical mean value of E@ as expected from 
Williams’ formula. On the other hand the experi- 
ments of Fowler and Oppenheimer*~ et al. on the 
scattering of electrons of mean energy between 6 
and 11 Mev in lead give only about 50 percent of 
the mean angle expected theoretically from 
Williams’ formula. Vargus has measured the 
scattering angle of some 55 particles of energy 
below 5X 10° ev and finds the mean value of E@ 
40 percent higher than the theoretical value. He 


1]. G. Wilson, Proc. Roy. Soc. A174, 73 (1940). 
*F. L. Code, Phys. Rev. 59, 229 (1941). 
+E. J. Williams, Proc. Roy. Soc. A169, 531 (1939). 

2 (sth Fowler and J. R. Oppenheimer, Phys. Rev. 54, 
5 N. L. Oleson, K. T. Chao, J. Halpern, and H. R. Cra 
Phys. Rev. 56, 482 (1939). eum = 
as 3) W. Sheppard and W. A. Fowler, Phys. Rev. 57, 273 


is not, however, sure whether this result is outside 
experimental error. 

In all these experiments, however, the method 
of measuring the energy of the particle is by 
measuring the curvature of its path in a magnetic 
field, and the scattering angle is measured by 
setting a cross-wire tangentially to the track of 
the particle below and above the scattering plate. 
Wilson! says that the accuracy of the measure- 
ments cannot be stated with certainty. In fact 
there is a systematic error due to particles (posi- 
tive or negative) being bent either in the same 
direction or in a direction opposite to the 
direction of scattering. The presence of a mag- 
netic field automatically excludes some low 
energy particles which are curved so much as not 
to reach the counter below the chamber. Besides, 
particles scattered through large angles are 
helped by the magnetic field to go out of the 
illuminated depth of the chamber. In fact, the 
magnetic field always sets a lower limit to the 
maximum scattering angle that can be observed 
in the chamber. 

Recently Shutt’ has measured the scattering of 
a large number of particles in 1 and 5 cm of lead, 
where he has not used a magnetic field. He does 
not know the energy of the particles photo- 
graphed, not even their energy distribution, very 


7R. P. Shutt, Phys. Rev. 61, 6 (1942). 
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accurately as he uses a 16-cm lead plate above the 
chamber, and this shifts the energy-spectrum by 
a considerable amount. Shutt has not tried to 
confirm Williams’ formula. He has assumed 
Williams’ equation and has calculated by an 
ingenious method the average nuclear scattering 
cross section of the mesons photographed. His 


value seems to be too low for a meson of spin 1, _ 


but agrees fairly well with the average value of 
Marshak and Weisskoff* if we assume Hartree’s 
energy distribution. In view of the theoretical 
calculations of Bhabha,® Weinberg,'® and Ma" all 
of which show a maximum of the cross section for 
a certain energy and then a falling off on both 
sides, I have tried to measure the scattering of 
mesons of different energy intervals in the way 
described below. The contents of this paper are 
concerned only with low energy mesons. 


EXPERIMENTAL 


The counter arrangement that controls the 
expansions of the chamber is shown in Fig. 1. 
Counter-pairs C,, C2, C3 are in coincidence and 
the three counters A (which are in parallel) are in 
anti-coincidence with C,C2C;. We can place lead 
of any thickness up to 6 cm between C; and A 
and also up to 25 cm between C, and C;. In the 
present experiment, the thickness of lead between 


C; and A is 6 cm, and there is no lead between C; . 


and C;. A 2-cm lead plate was used as the scat- 
tering material in the first part of the experiment 
and a plate 4 cm thick in the second part. The 
9-cm lead plate placed above the top counters is 
used only to cut out the soft component. From 
the arrangement it is obvious that mesons of 
range more than 2 cm but less. than 8 cm are 
photographed in the first series, and those having 
a range between 4 and 10 cm are photographed 
in the second series. The lead plate above the 
chamber excludes electrons and slows down some 
high energy mesons. From the range energy rela- 
tion of mesons as given by Rossi, the kinetic 
energies of the mesons photographed in the first 


® R. E. Marshall and V. K. Weisskoff, Phys. Rev. 59, 130 


1941). 
: 9H. J. Bhabha, Proc. Roy. Soc. A178, 314 (1941); Phys. 
Rev. 59, 100 (1941). ° 

10 i: W. Weinberg, Phys. Rev. 59, 776 (1941). 

uS,. T. Ma and C, F. Hsuh, Proc. Camb. Phil. Soc. 40, 


171 (1944). 
2B. Rossi and K. Greisen, Rev. Mod. Phys. 13, 249 


(1941). 


series lie between 0.55 X 10* ev and 1.3X 108 ev, 
and one meson of the second series lies between 
0.8 X 108 ev and 1.55 X 10° ev. The energy interval] 
in the two cases does not differ much. By placing 
larger amounts of lead between C; and C3, we can 
make the energy interval shift to higher values, 
That is, if the thickness of lead between C, and C, 
is x cm, the particles photographed will have an 
energy corresponding to a range greater than x 
but less than x+6 cm. It is the intention to use 
this method to measure the scattering of mesons 
of higher energies also. 

Before starting the experiment, the rate of 
counts C,C2C; without the anti-counters was first 
recorded. This gives the rate of mesotrons of all 
energies striking the scatterer. Next the efficiency 
of the anti-counters was tested in the following 
way. No lead was placed between the counters C, 


Fic. 1. Counter arrangement. 


Le, 


TABLE I. Counting rates. 


Thickness of 
Counting Rate of counts lead between Thickness of 
arrangement per hour CsandA scatterer 
28.2+0.75 2cm 
0.25+0.03 0cm 2cm 
1.45+0.07 6 cm 2cm 
1.58+0.08 6cm 4cm 


and A and the rate of counts with the anti- 
counters on (i.e., CiC2.C;—A) was found to be 
about one in four hours which is less than one 
percent of the total rate. The 6-cm lead plate was 
now placed between C; and A and the rate of 
coincidence C,;C2C;—A was noted for series one. 
This rate comes to 5.1 percent of the total rate. 
The rate slightly increases in the second series 
when the scatterer inside the chamber is 4 cm 
(see Table !). 

' It will be seen that since we do not use a 
magnetic field it is easier to measure the angle 
between the tracks of the particles above and 
below the scattering material, and at the same 
time we know that the energy of the particles lies 
within certain limits. The scattering angle was 
measured by fixing each negative between two 


‘clean glass plates on the stage of a very low 


power (magnification 2) microscope. This stage 
could be rotated. A small mirror was fixed to the 
periphery of this stage, and the angle of rotation 
of this stage was measured by a telescope and 
scale arrangement. As the position of the tracks 


‘in the negatives was different from different 


photographs, the two glass plates holding the 
negative between them were bodily moved until 
the point of intersection of the track above and 
below the lead plate was on the axis of the 
microscope. At this position the negative was 
rigidly fixed, and then by rotating the stage, a 
very fine cross-wire (which was sufficiently long 
to cover the whole negative) in the eyepiece was 
made to coincide first with the track above and 
then with the track below the scatterer. Half the 
deflection obtained in the scale divided by the 
distance of the scale from the mirror gave the 
actual angle of scattering. The accuracy with 
which each setting could be done was about 10’, 
and the mean of several settings should be 
accurate to 0.1°. 

The defect of the above method lies principally 
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in the fact that we do not know accurately the 
energy of each particle, but only know that it lies 
in a certain energy interval. Hence the compari- 
son with the theory could only be made in a way 
slightly different from Wilson and Code.’ 


COMPARISON WITH THEORY 


The theoretical development of electrical scat- 
tering of mesons and electrons is mainly owing to 
Williams.* He has taken into account screening 
of extra-nuclear electrons, and the arithmetic 
mean deflection for a finite nucleus is given by 


J=(19,.5—3.1 logio 2) 48, (1) 
(2) 


Mc*é represents the total energy of the incident 
particle, for = 1/(1—?)!. If E denotes the total 
energy, then dividing both sides of (1) by mc’, we 
get upon substituting for 6 


J  (19.5—3.1 logie Z)*(Nt)*Ze?* 
2mc? Ep? - mc? 


Substituting for e?/mc? (classical radius of the 
electron) and for N (number of atoms per cc) and 
Z, we have 


EB? = 0.906(4/t/@) X 10° ev, 


where @ is the average deflection in degrees. 

Now the mean kinetic energy of the mesons 
that we have photographed is 108 ev, and hence 
the mean total energy is 2 10°. The mean value 
of 6? that we can take is, therefore, given by 


2-10 108/(1—6*)!, 
B?=3/4. 
Substituting for 6? we have finally 
E=1.208(+/t/6) X 10° ev. (3) 


We shall now see how our results fit in with this 
theory. All the photographs were carefully 
analyzed, and only those tracks which have been 
recorded caused by the particles being stopped in 
the lead plate below were taken for measurement. 
Particles which made such large angles with the 
vertical that they were outside the solid angle 
covered by the anti-counters were excluded from 
measurement. Angles of scattering of individual 
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Fic. 2. Scattering in lead. 


No. of particle scattered in 2cmPb 


tracks being measured were classified into groups 
of 0° to 1°, 1° to 2°, 2° to 3°, and so on. Figure 2 
shows the distribution of particles with scattering 
angle for t=2 cm, and Fig. 3 shows the same 
for cm. 

The mean deflection for ¢=2 cm is 4.14°+.08, 

for ¢=4 cm is 5.67°+.10. 
Equation (3) shows that if E remains constant, 
which is very nearly true in our case, 0//t should 
be constant for the two thicknesses. The values 
of 6/s/t for t=2 and 4 cm are 2.93 and 2.83, 
respectively. The difference in these values is 
about 3 percent only for an increase in +/t by a 
factor 1.4. Considering that for t= 4, the value of 
_ Eisa little greater than for ¢=2, we can say that 
the +/t law is well confirmed. 

Let us now see how the-mean deflection agrees 
with the energy. If we put @=4.14° and ¢=2 cm 
in Eq. (3), we get the value of E to be 4.13 X 10’ ev. 

The maximum value of £ possible for the 2 cm 
series is, however, only 2.3 X 108 ev, the minimum 
being 1.55108 ev. If we take the energy spec- 
trum in this interval to be smooth, the mean 
value of the total energy E of the particles 
registered is actually 1.93108 ev which is only 
47 percent of the value found out from William’s 
formula. Consequently we may say that the 
average deflection given by Williams equations is 
53 percent too large. This is not in agreement 
with Wilson! or Code,? who have verified Williams’ 
equation for high energy mesons. Fowler and 
Oppenheimer*-* al., have, however, found 
almost the same amount of deviation for electrons 
of mean energy between 6 and 11 Mev, as we 
_ have found for mesons of low energy. Since the 
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Fic. 3. Scattering in lead. 


mean angle of multiple scattering for high energy 
mesons satisfies Williams’ equation whereas the 
low energy mesons show a marked deviation 
from the theory, the cause of the discrepancy 


seems to depend on the energy of the mesons - 


concerned. We have seen that the 7/¢ law is 
confirmed satisfactorily, and hence the statistical 
effect of all the nuclei taken together agree with 
the theory. It therefore appears that the electrical 
scattering of each individual nucleus is to be 
modified in such a way as to give a smaller mean 
angle for lower energies than the theory. 


NUCLEAR SCATTERING 


We can see from Fig. 2 and Fig. 3 that the 
distribution of particles with angle of scattering 
shows a greater amount of deviation from 
Gaussian distribution for t=4 cm, than for t=2 
cm. This is to be expected since the nuclear 
scattering cross section is proportional to the 
thickness of the scatterer, whereas multiple scat- 
tering varies as \/t. Undoubtedly, the percentage 
of single scattering, which is generally large angle 
scattering, is more in Fig. 3 than in Fig. 2. Shutt 
has developed a method of calculating the aver- 
age nuclear scattering cross section from the 
scattering data for two thicknesses of lead. The 
main idea is the following. If J(@, t)d@ represents 


the multiple scattering intensity between @ and 


6+d0, expressed as the percent of the total 
number of traversals, then instead of calculating 
the distribution with respect to @ we transform @ 
to a new variable u given by 


u=6//t.. 
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We then have 
f(u) =I(8, t)(00/du) = I(8, t)/t, 


where f(u)du is the number of particles in the 
range between u and u+du. If the scattering is 
entirely multiple, the f(u) values for all angles 
would be identical for the two thicknesses of lead 
except, however, for statistical fluctuations. 
Following Shutt we have calculated the f(u) 
numbers, in the range of u from 2 to 12 for both 
the thicknesses, and they are given in Table II. 
Summing up the f(u) numbers we find that in- 
stead of an equality the value for the 4 cm thick- 
ness is greater than that for 2 cm by 


6=116.8—56.4=50.4. 


The method used by Shutt for calculating the 
number of singly scattered particles is not strictly 
correct. The difference between the figures in the 
second and third columns of Table II shows that 
multiple scattering due to the electric charge of 
the nucleus is not the only scattering that pre- 
vails. We therefore proceed as follows. We assume 
that J(@,¢) represents the sum of multiply and 
singly scattered particles J,, and J,, respectively, 
so that 
I(0, t) =In+I,. 


I, is directly proportional to the thickness of the 
scatterer, and consequently we can write 
I, = kt, 


where k gives the percentage singly scattered in 
1 cm of lead.” That is 


I(0, t) =In+ht, 
or 


f(u) =1(8, t)/t= (Im +Kt)V/t = In(v/t) 


Since for multiple scattering J,,./t remains con- 
stant for different values of t, 6 the difference 
between the values of f(u) for ‘=4 and ¢=2 will 
be given by 

6=K(4!—2!) =5.172K. 


The experimental value of 6 is 50.4, which gives 
K=9.73+2.4 percent. 


The statistical error has been calculated in the 


13 is a function of @ but not of ¢t. The K that occurs in 
- ce for 6 further down is the integral of this & over 
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TABLE II. Values of f(u). 


degrees x cm? % x %x 
1 
5 
8 
6 
4 


same way as Shutt. According to Shutt a correc- 
tion of 30 percent is to be added to this for using 
the projections of the angles in the plain of the 
cloud chamber instead of the actual angles in 
three dimensions. Making this correction, the 
single scattering in one centimeter of lead comes 
to 
K=12.65+3.2 percent. 


The average cross section ¢ per nucleon is given 
by 

o=K/N, 
where N is the number of nucleons per cubic 
centimeter of the material, since K is the per- 
centage scattered in 1 cm of the material. 
Inserting values we get 


12.65 X 10-* 


= 10-* 
6.024 X 11.4 


+25 percent per nucleon. 


This is about twenty-five times the average value 
obtained by Shutt for the whole energy distri- 
bution. Bhabha® has calculated the nuclear scat- 
tering cross section for tranversely polarized 
mesons taking radiation reaction into account 
according to the classical theory. Taking the con- 
stant g2*x?/h = 1/13.3, as is required by the theory 
of nuclear forces, he finds that the cross section 
attains a maximum of 3X10-** for E=3.5yuc? or 
3.5108 ev. The mean value of our small energy 
spectrum is 2X 10* ev, and the value of the cross 
section found in this experiment agrees very well 
with Bhabha’s theoretical value 1.6 10-** for 
this energy. Recently Ma and Hsueh" have 
calculated the nuclear scattering cross section for 
transverse mesons according to the quantum 
theory, and their value is 2X10-** cm?* per 
nucleon for E=2X10* ev. Weinberg’s'® theo- 
retical value of the cross section for mesons of 
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Fic. 4. Track of ionizing particle scattered in a 4-cm 
lead plate. 


spin } is of the same order of magnitude as ours, 


and agrees better with our value than with 
Shutt’s. The experimental value of the cross 
section obtained by Shutt should be necessarily 
lower than ours since his value is an average for 
the whole energy spectrum. The presence of high 
energy mesons in large numbers completely 
masks the higher value of the cross section for 
low energy mesons. This is the reason why 
Wilson! ‘or Code? (as calculated by Shutt’), who 
have deliberately excluded mesons of energy 
below 2X10* ev, get values for the nuclear 
scattering cross section which are too low. Here 
on the contrary, we have excluded the high 


energy mesons, and hence the value that we 
obtain agrees better with theoretical values for 
low energy mesons. 


SCATTERING OF A PROTON 


Figure 4 shows the track of a heavily ionizing 
particle being scattered through an angle of 5.6°, 
in a 4-cm lead plate. The track cannot be that of 
a meson, for if it were a meson slow enough to 
produce such large ionization it would not have 
penetrated the 4-cm plate of lead. 

The range of the particle lies between 4 and 10 
cm of lead and hence the minimum kinetic energy 
possible for the particle, taking it to be a proton 
is 10° ev, or the total energy is 10° ev. The mean 
angle of multiple scattering for this energy, ac- 
cording to Eq. (3) is 2.4°. The probability that a 
proton of this energy be scattered through 5.6° 
(the observed angle) is 7.3 percent. It will be still 
lower if the range is nearer to 10 cm. So the 
scattering seems to be more probably a nuclear 
one. There is another alternative. The charge of 
the particle may be more than one. If the inci- 
dent particle is of charge Ze, the right-hand side 
of Eq. (3) is to be multiplied by Z. If we take the 
particle to be doubly-charged, the mean angle 
would come to 4.8° which is not very different 
from the observed angle. This alternative cannot 
be excluded as the ionization produced by the 
particle is almost that of an a-particle and this is 
what would be expected for a doubly-charged 
proton. 

The author wishes to express his grateful 
thanks to Professor H. J. Bhabha for suggesting 
the experiment and his inspiring guidance 
throughout the work. He also wishes to thank 
Professor Sir,C. V. Raman for lending the appara- 
tus used for measuring the scattering angles. 
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The interaction of slow or C neutrons can be modified by changes in the physical parameters 
of the material which the neutrons are traversing. Measurements of the effect of grain size on 


total C neutron cross section were made on randomly-oriented poly-crystalline iron and copper. 
In the case of copper it was found that the total cross section decreased in an exponential 
manner with increase in grain size. For iron it was simply established that the total cross 


section was substantially smaller for the larger grain size. 


HE work of Nix, Beyer, and Dunning! has 

shown that the presence of order in binary 
iron-nickel alloys, in- which the grain size was 
maintained constant and free of preferred orien- 
tation, decreased the cross section below that 
obtained from random alloys: the greater the 
degree of order, the lower the total cross section 
obtained. The establishment of such an effect"? 
illustrates that thermal neutrons may be used 
as a tool to study the structure of the solid state. 
In order to render this tool more useful for our 
systematic studies of order-disorder transforma- 
tions in alloys, it was necessary to establish the 
influence of other physical parameters on the 
neutron scattering. As our first problem in this 
series of studies, we selected the effect of grain 
size. Grain growth frequently takes place in the 
thermal treatments necessary to produce order- 
disorder transformation. As metals we selected 
copper as a non-ferromagnetic face-centered 
cubic metal free from polymorphic transforma- 
tions, and iron as a ferromagnetic body-centered 
cubic metal subject to polymorphic transforma- 
tions. 


EXPERIMENTAL ARRANGEMENT AND.- 
PROCEDURE 


The neutron source was a radium-beryllium 
bulb containing 200.5 mg of radium and 8.0 mg 
of beryllium. Figure 1 shows the detailed arrange- 
ment of the “howitzer,”’ scatterer, ionization 
chamber, amplifier, and counting circuit. Ther- 
mal neutrons from the “‘howitzer’’ passed down 
a cylindrical collimating tube 10 cm in diameter. 

* Now with the Research Laboratories of The Sharples 
pate. Philadelphia, Pennsylvania. 

1F. C. Nix, H. G. Beyer, and J. R. Dunning, Phys. Rev. 
58, 1031 (1940). 


(1940) G. Beyer and M. D. Whitaker, Phys. Rev. 57, 976 
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The interior wall of the tube was lined with 
cadmium. The collimating tube contained a 
0.635-cm hollow outer shell which was filled 
with B,C in order to reduce the number of 
neutrons escaping into the room. The length of 
collimating tube from source to test sample was 
27 cm. The collimating tube was constructed in 
two parts permitting samples to be placed in the 
neutron beam at the center of the system. The 
collimating tube between sample and detector 
was 26 cm long which prevented the C neutrons 
scattered at an angle of greater than 8 degrees 
from reaching the detector. 

The neutron detector was a BF; ionization 
chamber connected to a linear amplifier and 
scaling system recorder located in a distant room 
in the laboratory building. 


Fic. 1. Experimental arrangement. 
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The samples to be tested were automatically 
introduced into and removed from the neutron 
beam by electric motors. The sequence of opera- 
tion of the motors was controlled by a timing 
circuit which providéd for the following neutron 
beam conditions: 

1. Open beam. 

2. Cadmium slab in beam. 

3. Sample and cadmium slab in beam. 
4. Sample in beam. 


When comparing the relative transmission of 
two samples we use the 1, 2, and 4 conditions 
with the cadmium replaced by the second sample. 
The 1 condition permitted us to keep a frequent 
check on amplifier and scaling circuit conditions. 
The timing circuit actuated the motors to oper- 
ate at intervals varying from 3 to 60 minutes. 
All measurements in this research were carried 
out at 80-minute cycles, which allowed 20 min- 
utes for each part of the cycle. 
The data were recorded automatically by 
photographing simultaneously the counter, an 
electric clock which could be read to 0.1 second, 
and signal lamps. The use of the clock permitted 
the duration of each part of the cycle to be 
accurately measured. The combination of the 
three signal lamps gave positive information as 
_to which of the above-listed beam conditions 
prevailed at any instant. A 16-mm motion pic- 
ture camera, provided with a feature allowing 
the exposure of one frame at a time, was used to 
photograph thé data, ‘which -included clock, 
signal lamps, and counter. Exposures were taken 
at intervals of 45 seconds. 
The thickness of the scattering sample was 
such that it removed approximately 40 percent 
of neutrons from the ‘‘open’”’ beam. 


PREPARATION OF SAMPLES 


The copper samples were prepared from an 
ingot of very pure deoxidized high conductivity 
copper by hot-working plates of the material at 
a temperature of 1000°C to 50 percent of their 
initial thickness of 3.05 cm, cold working to 
a thickness of 0.76 cm, and milling to a final 


thickness of 0.63 cm. The samples were plates 


12.7 X 12.7 X0.63 cm. 

The various grain sizes were obtained by an- 
nealing the cold worked material. The annealed 
iron sample was obtained by vacuum annealing 
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at 1000°C for 30 minutes. Seven copper samples 
were prepared by vacuum annealing for 30 min- 
utes at each of the following temperatures, 
450°C, 500°C, 550°C, 650°C, 750°C, 800°C, 
850°C, while an eighth sample was annealed at 
1000°C for one hour. Cold worked samples of 
iron and copper were also retained. 

_ The grain sizes of all scattering samples were 
obtained from photomicrographs of the polished 
and etched specimen. The qualitative orientation 
of the crystallites in each sample was established 
from x-ray diffraction patterns taken at grazing 
angles. The x-ray and microscopic examinations 
were made on the portion of the specimen ex- 
posed to the neutron beam. 


METHOD OF MEASUREMENTS 


a. Relative Neutron Transmission 
Measurements 


The relative neutron transmission measure- 
ments simply consist in obtaining the transmis- 
sion of two identically-shaped plates in the 
neutron beam. The number of counts per minute, 
Na, registered in the ionization chamber with 
annealed sample as scatterer, is compared with 
the number obtained with cold worked sample 
as scatterer, V5. 

The transmission of the annealed sample rela- 
tive to the cold worked sample, 74», is expressed 
by, 

Ta= X 100. (1) 


The final value of relative transmission was 
taken as the average of 36 cycles. 


b. Measurement of Total Cross Section 
of Scatterer 


The total cross section is obtained from the 
following four measurements: 


1. The number of disintegrations produced in the BF; 
ionization chamber per minute, Ni, by the open beam 
with no scatterer present. 

. The number of disintegrations per minute, Ne, with 
0.45 g/cm? of cadmium in the neutron beam. 

. The number of disintegrations per minute, Ns, with 
both cadmium and the scatterer whose cross section is 
to be determined in the beam. 

. The number of disintegrations per minute, N,, with the 
scatterer alone in the neutron beam. 


was 


The cadmium in parts 2 and 3 of the cycles 
serves to prevent any thermal neutrons from 
reaching the ionization chamber. 

The transmission, 7, of the scatterer for C 
neutrons can then be determined from the follow- 


ing equation, 
T= (2) 


The total cross section of the scatterer can 
be calculated from the equation, 


(3) 


where a is the total cross section in sq. cm; m, the 
number of nuclei per cubic cm of scatterer; and 
t, the thickness of the scatterer in cm. The 
number of nuclei per cm* can be obtained from 


the formula, 
n=dX6.02 X 10*8/Atomic weight, (4) 


where 6.02 X 10** is Avogadro’s number and d is 
the density of the scatterer. 

The cross-section values given in Tables I and 
II were taken from the average of 18 determina- 
tions. The number of neutrons recorded in the 
ionization chamber in the open beam condition 
was approximately 1000 per minute. 


DISCUSSION OF RESULTS 


The percent changes in transmission, referred 
to the cold worked sample, for measurements on 
iron and copper samples are given in Table I. 
These values are obtained from Eq. (1). They 
show the annealed samples to possess the larger 
transmission ; with random orientation of crystal- 
lites, the larger the grain size the greater the 
transmission. 

The total cross-section values, ¢, obtained on 
cold worked and annealed samples are given in 
Table I]. These values were obtained by the 


TaBLE I. Relative thermal neutron cross sections. 


Grain Percent change 

size in in transmission 

Preparation milli- Orienta- _ referred to cold 

Material of sample meters tion worked sample 


Fe Cold worked 0.077 +7+42.0 
Fe Annealed$hr. 0.127 Random 

1000°C 
Cu Cold worked 0.005 +9+2.0 
Cu $hr. 0.081 Random 
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TABLE Il. Thermal neutron cross section as a function of 
grain size X 


Grain Total 
size in cross-section ¢, 
tion milli- Orienta- x<10™ 
Material of sample meters tion cm 


Fe Cold worked 0.077 Random  13.06+0.4 
Fe Annealed$hr. 0.127 Random  11.39+0.3 


Cu Cold worked 0.005 Random  10.84+0.3 
Cu Annealed$hr. 0.017 Random  10.16+0.4 


Cu Annealed$hr. 0.015 Random  10.13+0.5 
500°C 


Cu 0.021 Random 9.96+0.4 
Cu $hr. 0.036 Random 9.62+0.3 
Cu $hr. 0.373 Random 7.76+0.2 


Cu Annealed$hr. 0.567 Random 7.58+0.3 
800°C 


Cu $hr. 0.640 Random 7.37404 
Cu Annealedihr. 2.40 Random 7.04+0.3 
1000°C 


use of the Eq. (3) from transmission values 
obtained from Eq. (2). In Table II is also in- 
cluded information on preparation of samples, 
values of grain sizes determined microscopically, 
and comments on the orientation of the crys- 
tallites. 

A plot of total C neutron cross section, ¢, 
versus grain size for copper is given in Fig. 2. 
The larger the grain size, the smaller the otal 
cross section. It is to be noted from Table II 
that all samples were free from any detectable 
preferred orientation of the crystallites. These 
results clearly show that the total C neutron 
cross sections of solid material cannot be looked 
upon as physical constants of the elements or 
compounds in question but are largely deter- 
mined by the diffraction effects within the body 
of the scatterer. With experimental arrangement 
given in Fig. 1 each Bragg reflection tends to 
remove neutrons from the well collimated beam, 
causing them to be deflected to the cadmium- 
lined collimator tube wall and consequently 
removing them from the beam. For completely 
randomly oriented crystallites, the greater the 
number of crystallites intercepting the beam the 
greater the decrease in the transparency of the 
sample. The data in Tables I and II and the 
graphical illustration of the copper data in Fig. 2 
clearly show the increase in transparency with 
the consequent decrease in total cross section 
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Fic. 2. Grain size vs. cross section X 10% cm~? for copper. 


with increasing grain size. With random orienta- 
tion, the larger the grain size, the fewer grains 
are available for Bragg reflection. 

These experiments clearly show that the 
thermal neutron cross-section values given in 
literature lose much of their supposed quanti- 
tative significance and should be regarded as 
purely qualitative in nature. As an illustration, 
the values for copper depend to a large extent 
on the grain size and would be further compli- 
cated by the presence of preferred orientation. 
In the case of alloys there is the additional com- 
plication of atomic distribution in solid solutions 
and spatial distribution of phases for hetero- 
genous binary alloys. In the work of Nix, Beyer, 


and Dunning! on iron-nickel alloys in the region 
around 75 atomic percent nickel, it was shown 
that the presence of ordering tends to increase 
the transparency to thermal neutrons. In this 
work the grain size was maintained constant 
while varying the degree of order, and there was 
demonstrated an increase of as much as 8 percent 
in transparency on increasing the degree of order. 

As stated in the introduction, this investiga- 
tion was undertaken with the object of exploring 
the possible use of neutrons as a tool for investi- 
gating the solid state. In order to evaluate 
properly data obtained at elevated temperatures 
on order-disorder or other solid-state transforma- 
tions, it was found to be necessary to establish 
the effect of parameters such as grain size and 
change in orientation of crystallites on neutron 
cross-section values. 

It appears that the use of thermal neutrons 
with wave-lengths comparable to the inter- 
atomic spacings in crystals will be a very useful 
method for studying the structure of the solid 
state. Accurate quantitative information per- 
mitting unequivocal interpretation of changes 
in the solid state must await the development 
of intense sources of monochromatic neutron 
beams; however, much valuable information can 
be obtained from the use of a relatively weak 
radium-beryllium source of wide wave-length 
bands such as were used in the present investi- 
gation. . 
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The binding energy of the deuteron and the scattering cross section of the proton by fast 


neutrons are calculated by using new forms of nuclear potential suggested by K. C. Wang. 
Results obtained are found to be in good agreement with experimental values when “zero 


cut-off” of the potential is employed. 


I. INTRODUCTION 


T has been pointed out by-K. C. Wang! that 
the force between two nuclear particles may 
be related to the gravitational force. He con- 
siders two alternative forms of the nuclear 
potential : 
V=—Ae*, (1a) 
and 
V=—(B/r)e*'", (1b) 
with 
k=h/mc=3.84X10- cm, 
where the constants A and B, as determined by 
the gravitational constants, are 4.78 x 10-“ and 
1.83 10-5, respectively. 

The purpose of the present work has been to 
determine whether the potential (1a) and (1b) 
can be used to obtain satisfactory values for the 
binding energy of the deuteron and the scattering 
cross section of the neutron by the proton. The 
calculations follow closely those of Bethe and 
Bacher.’ 


Il. THE BINDING ENERGY OF THE DEUTERON 


The wave equation for the relative motion of 
the two nuclear particles is 


Ay+(M/h*)(E—V)y=0. (2) 


The potential is spherically symmetric. Equation 
(2) can thus be separated in polar co-ordinates 
r, 6, , by putting 


V(r, 8, 6) 
where Pim is a spherical harmonic. The wave 


1K. C. Wang and H. L. Tsao, Phys. Rev. 66, 155 (1944); 
Nature i55, 512 (1945). 
936). Bethe and R. E. Bacher, Rev. Mod. Phys. 8, 82 
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equation for u; is then 


For the ground state /=0, we have 
V)uo=0. (3) 


To solve (3), two forms of potential cut-off given 
by (1a) and (1b) are employed; 


(a) zero cut-off: V=0 for r<a, V= V(r) for 
r>a; 

(6) straight cut-off: V=V(a) for r<a, 
V= V(r) for r>a. 


(a) Zero cut-of:.—For r<a, V=0; Eq. (3) 
takes the form, 


(d*uo/dr*) —(Me/h?)uo=0, 


where e= — E, the binding energy of the deuteron. 
Its solution is 


uo =D — } 
where 
B=(Me/h*)!, 
and 
j duy/dr = DB(e*" +-e-*"). 
Therefore, 
1 duo 
(— = 8 coth Ba. (4) 
Uo dr 


For r>a, V= V(r); the asymptotic solution of 
Eq. (3) is 
Uo =ce~Fr 4 
Treating c as a slowly varying quantity whose 
second derivative with respect to r may be set 


gion 
ease 
this 
tant 
was | 
cent 
tiga- 
ring 
esti- 
uate 
‘ures 
blish 
and 
tron 
rons 
iter- 
seful 
per- 
nges 
nent 
i 
can 
veak 
ngth 
esti- 
4 
4 
|_| 


164 MU-HSIEN WANG 


we hove For r>a, [(1/o)(duo/dr) is the same ag 
duo dc duo given by (5). Therefore, instead of (6), we have 

dr dr dr? cot ya. (6a) 


Substituting this last expression in (3), we get the 
Equations (6) and (6a) serve to determine the 


Taste I. Values for Eq. (7). binding energy of the deuteron, if we use the 
experimental value of — V(a) = Vo (the depth of 
cath the potential well) and calculate a (the range of 
3.50 1.33 1.38 the nuclear force) from (1a) or (1b). Conversely, 

3.60 1.32 1.32 
3.70 131 1.26 if we use the experimental value of ¢, we can 
calculate the value of Vo and a by the aid of the 

Eqs. (1a) or (1b). 

differential equation for c, ' (c) Numerical calculation.—(1) V=0, for r <a; 
: V= V(r) = —Ae*'", for r>a. We take the recent 
dc experimental value, Vo= 10.5 Mev, then a calcu- 
dr 26h? lated from (1a) is equal to 4.2110-* cm, 


Equation (6) becomes numerically 


the solution of which is 
coth [1.24 10"*(Vo—e)#a]=(Vo/2e)—1. (7) 


c= Fexp (- Soe r). From Table I we see that the root of (7), the 

Thus binding energy of deuteron, is 3.60 10~- erg or 
2.26 Mev. i 

F exp (2) V=0, for r<a; V(r) = —(B/r)e"*, for r>a. 

and Also take Vo= 10.5 Mev, then a=4.42 X 10-* cm 

1 duo _M M|V@| V@|_ from (1b). 
ye . (S From Table II we get the binding energy of the 
deuteron, 3.66 X 10-6 erg or 2.30 Mev. The results 


In order that the wave function uo could be given by (1) and (2) are in good agreement with 


joined smoothly at r=a the expressions (4) and _ the experimental value 2.27 Mev. 
(3) V=V(a), for r<a; V(r) = forr>a. 


(5) must be equal; i.e., 
Also take 10.5 Mev, then a=4.20X 10-" cm 
M|V(a)| from (1a). Equation (6a) becomes numericall 
—_ ©) 0.5Vo/9/ 
cot 
Straight cut-off.—For r<a, V=V(a); Eq. (Vo—«)3 
(S) takes the form, ' TABLE II. Values for Eq. (7). 
Vo—e)uo=0. —coth  Ve/2e—1 
3.50 1.29 1.38 
— 3.60 1.29- 1.32 
Its solution is 3.70 1.27 1.26 
3.80 1.27- 1.19 
with 0 Y’; 
(4) V=V(a), for r<a; V(r) = —(B/r)e*", for 
Fa (1V(@)| r>a. We get the results listed in Table III. 
Thus From Tables III and IV, the binding energy of 
(— 1 duo inincitiicein the deuteron is 8.0X10-* and 8.3X10-* erg, 
uo ar —* respectively. These values are too large in com- 


era 


com- 


TABLE III. Values for Eq. (6a). 


0.5Vo, 
«,10erg cot [1.24 ( 
7.8 0.018 0.067 
8.0 0.038 0.037 
8.2 0.061 0.027 
TABLE IV. Values for Eq. (6a). 

«, 10-* erg cot [1.24 ( 
8.2 —0.020 0.206 
8.3 —0.013 0.005 
8.4 —0.004 —0.083 


parison with the experimental value. The correct 
results will be obtained if we take Vo=5.93 Mev. 


Ill. PROTON-NEUTRON SCATTERING 


Let us denote by E the kinetic energy of a 
proton and a neutron in a coordinate system in 
which the center of gravity of the particles is at 
rest ; this is equal to one-half of the kinetic energy 
of the incident neutron in a system at rest. The 
radial function 1; will satisfy the equation, 


2/d2u, I(l+1 


Asymptotically for large r, the solution of (8) is 
sin (kr—3lx+5,) (9) 
k?= ME/h’. 


with 


Then the cross section de i.e., the number of 
neutrons scattered per unit time through an 
angle between @ and 6+4d4, if there is one neutron 
crossing unit area per unit time in the incident 
beam,—is given by the well-known formula 
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It has been shown? that if (1/k)>a, a being the 
range of the force, all phases 6; will be small 
except 59. Then 


do =2xk-? sin? 59 sin 
and 


o= sin? do. (11) 


For the ground state of the deuteron, we have 
already shown that 


dug~ MV, 


Now,’ in the present case E is positive, we should 
have*® 


1 duy M(E +e) 
) uotug dr, 
dr r=a h?ugtug- 0 
M(E+¢ 


where and are uo functions for r<a, 
corresponding to the positive and negative values 
of E respectively ; i.e., 


ot=Bsinkr, =D{e"—e-*}, 


thus 
e 8 coth Ba—k cot ka 
(14) 
For large r, 


1 dugt 
+ dr 


) =k cot (ka+ (15) 


In order to join the wave function smoothly at 
r=a, the expressions given by (13) and (15) 


‘must be equal, thus 


k cot =A. (16) 
3 Cf. reference 2, page 115. 


. 


TABLE V. Numerical results. 


Ei Mev ¢ ¢ (obs.) 
2.15 0.85 10-** 1.43 x 10-*4 1.2 10-*4 10-*4 0.5-0.8) 10-** 
1.05 1.03 x 10-*4 1.86 x 10-*4 10-*4 2.4X 10-** 1.1-1.5)10-** 
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Under the assumption that ka is small, we get 
cot d9=A/k, (17) 
and by (11), 
o =4n/(k?+A?). (18) 
If 9 is small, A ~a, we have 
oo (19) 


The numerical results are given in Table V, in 
which we take V = — Ae*’", and the values of a are 
obtained from the Table I at the point of joint. 
Corresponding to (18) and (19), Bethe and Bacher 
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have derived the formulae: oog=4mh?/(e+E), 
op and the values calculated from these 
are given in the 4th and 5th columns of Table V 
for comparison. 

It is seen that the present results are in better 
agreement with the experimental results than 
those given by Bethe and Bacher. 

If we take V= —(B/r)e*/", the results do not 
differ appreciably from those given above. 

In conclusion, the author wishes to express his 
thanks to Dr. K. C. Wang for suggesting this 
calculation and for helpful discussions. 
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field and of the wave fields of elementary particles are obtained with the help of the invariant 

functions of Jordan, Pauli, and Dirac, together with their generalizations. It is shown that the 

Fourier transforms of these potentials are closely related to those of the outward and inward 

moving waves given by Dirac for the scattering problems in quantum mechanics, and their 

connection is discussed. It is also shown that there exists a type of potentials which represents 
_ waves with frequencies of opposite signs propagating in opposite directions. 


1. INTRODUCTION 


NUMBERS 7 AND 8 OCTOBER 1 AND 15, 1945 
Fourier Transforms of Retarded and Advanced Potentials 
S. T. Ma 
National University of Peking, Kunming, China 
(Received March 20, 1945) 
The Fourier transforms of the retarded and advanced potentials of the electromagnetic 
with 
=—-V?. 3 
(3) 


ORDAN and Pauli’ introduced into the quan- 
tum theory of electromagnetic field a function 


5(xo+ |x|) —8(xo— |x|) 
|x| 


which plays a very important part in the commu- 
tation relations of the field variables. Following 
Dirac,? we shall denote it by —A and write © 


5(x0— |x|) —5(xo+|x|) 

|x| 
where X denotes the four-vector (x, xo). This A 
function is relativistically invariant and satisfies 


the wave equation of electromagnetic field in 
vacuum 
(2) 


A(X) = (1) 


OA(X) =0, 


! Jordan and Pauli, Zeits. f. Physik 47, 151 (1928). 
? Dirac, Proc. Roy. Soc. A180, 1 (1942). 


It can be expressed as a triple Fourier integral of 
the form 


where k= |k|, or as a quadruple Fourier integral 
of the form’ 


A(X) f A(K) exp (—ikero) 
Xexp (+7k-x)dtk, (5) 


where K denotes the four-vector (k, ko). 

Some functions of the same type but more 
general than that of Jordan and Pauli appeared 
first in Dirac’s‘ positron theory, which are in the 


3 Dirac, Ann. I'Inst. H. Poincaré 9, 13 (1939). 
* Dirac, Proc. Camb. Phil. Soc. 30, 150 (1934). 


1945 


(3) 


(4) 


ore 


the 


notation of recent papers by Pauli,’ 


1 sin (k?-+ «2x 
Xexp (+7k-x)d*k, (6) 
1 cos (k?+ x?x9)! 
D,(X)= 


(2x)* (k?-+«?)! 
Xexp (+7ik-x)d*k. (7) 


These functions are also relativistically invariant 
and satisfy the wave equations 


(0 +x?)D(X) =0, (8) 
(0 +x?)Di(X) =0. (9) 
As shown by Dirac, if we denote |x| by x, 
(10) 
( aie (x, Xo ’ 
where 
J —x?)*] 
F(x, Xo) =40 (x>xo> (11) 
— Jol «(xo?—x*)#] (xo<x), 
and 
D,(X) (x, Xo), (12) 
4m x Ox 
where 


(x>x, 
A(z, (x9>x, *) 3) 


(x>xo>—x). 

The function D(X) reduces to A(X)/4x when 

x=0 and may be regarded as the analog of 

4/4 in field theories of elementary particles 

whose rest mass does not vanish. The function 
D,(X) reduces when x=0 to 


1 
=— 


29? x? 


for Xo. (14) 


The possibility of employing the D, function in 
quantized field theories has recently been dis- 
cussed by Pauli. 

It is of interest to point out that besides the 
positron theory and the general theory of second 
quantization there isa third type of mathematical 
problem in which these D functions can be used 
with advantage; namely, the Fourier resolution 
of retarded and advanced potentials of field 
theories which will be carried out in the present 


5 W. Pauli, Phys. Rev. 58, 716 (1940); Rev. Mod. Phys. 
13, 203 (1941). 
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paper for both the electromagnetic field and the 
wave fields of elementary particles with the help 
of these D functions and their generalizations. 
The mathematical treatment will be developed 
on the lines of Dirac’s® treatment of outward and 
inward moving waves in scattering problems of 
quantum mechanics. The connection between the 
Fourier transforms of the retarded and advanced 
potentials and those of the outward and inward 
moving waves will be discussed. 


2. FOURIER TRANSFORMS OF RETARDED 
AND ADVANCED POTENTIALS OF 
ELECTROMAGNETIC FIELD 


Let A(X) be a vector function with four com- 
ponents representing the scalar and vector po- 
tentials of the electromagnetic field when the 
charge and current densities are represented by a 
vector function j(X). It is known in classical 
electromagnetic theory that the function A 
satisfies the inhomogeneous wave equation 


=47j(X) (15) 


and that there are two independent solutions of 
(15) which have simple physical meanings, 
namely, the retarded potentials 

j(x', 


A(X) = a°x (16) 


and the advanced potentials 


x— x’ 

|x—x’| 

These retarded and advanced potentials may be 

expressed as four-dimensional integrals with the 

help of Dirac’s 6 functions, thus 


Acx)= (18) 
or 
fice, a9) 


if we write r=x—x’, r= |r|. Now the 
6 functions satisfy the relation 


5(ro—r) + +r) = 8(ro?—r?) (20) 


2r 


* Dirac, Zeits. f. Physik 44, 585 (1927); Principles of 
Quantum Mechanics (Clarendon Press, 0 
second edition. 
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and (1) and (20) give us 


5(roFr) 


= 3A(R), (21) 


where R denotes the four-vector (r, 76). Formulae 
(19) may therefore be written 


A(X) = f Jax’. (22) 


It has been shown by Dirac that, if we write 


OrA(R) =0. 
Using Dirac’s method.we can also show that 
O1r6(ro?—r*) =475,(R), 
54(R) = 5(r2)5(r3). 
From (21), (23), (24) it follows that 


O 475,(R). 


(23) 


(24) 
where 


(25) 


The relations (23)-(25) correspond to the fact 


Now 


S. T. MA 


that the retarded potentials, advanced potentials, 
and their mean values satisfy the same inhomo. - 
geneous wave equation, while the differences of 
the retarded and advanced potentials satisfy the 
homogeneous wave equation for the vacuum and 
are free of singularity everywhere. 

Let the vector functions A(X) and j(X) be 


. expressed as Fourier integrals in the form 


1 
(20)? 


1 
(2m)? 


f A(K) exp (i(K, X))d*, (26) 


j(X) = 


where 


f j(K) exp (i(K, X))d*k, (27) 


(K, X) =Roxo—k-x. 
From (22), (26), and (27) we find 


{8(ro?—r?) +3A(R) } 
Xexp (—i(K, R))d‘r. 


Thus the Fourier transforms of the potentials can 
be expressed in terms of those of the source 
densities if we know the Fourier transforms of the 
functions 6(ro?—r?) and A(R). 


A(K) =3(K) f 
(28) 


f exp (—7(K, R))d*r= f f exp (—ikoro+ik-r)d*rdro 


2r 


so that, on writing 


ko 
cos kor exp ( ) a8, 


cos kor exp (k-r) | 
-f d*r, 
r 


(29) 


1 
4:(K) 


we have 


f exp (—i(K, 


r 


(30) 


From (5) applied with K and X replaced by R and K, we obtain 


} f A(R) exp (—i(K, R)*)dr = —2mi?A(K). 


Substituting from (30), (31) in (28), we obtain 


A(K) = 21*j(K) {4.(K)¥iA(K)}, 


(31) 


(32) 


a relation which expresses the Fourier transforms of the potentials in terms of those of the source 
densities with the help of the D functions. We have written here A, A, for the values of 41D, 4xD,, 


respectively, in the case of x=0. 


d*r, 


29) 


30) 


31) 


32) 


rce 
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3. A MODIFIED FORM OF RETARDED AND ADVANCED POTENTIALS 


Substituting from (1) and (14) in (32) we obtain for the Fourier transforms of retarded and 
advanced potentials, the expressions 


k) —5(ko+k) 
uk k 


The expressions on the right-hand side of a are not defined for ko= +k, where the term 1/(k?— ko?) 
tends to infinity. In order to give this term a definite meaning for ko = +k, we can replace the function 
1/(k?—ko*) by the function 7A,/2, defined by (29), at the points of singularity. When thus interpreted, 
(33) substituted in (26) gives our original formulae (19). 

There is, however, another way of dealing with the singularities at kj = +k. Following Dirac’s 
works on scattering, we can take an integral whose integrand contains 1/(k®—k»?) as a factor to be 
equal to its Cauchy principal value, e.g., 


exp ((K,X)) j(K) exp (i(K, X)) 


It will be interesting to see what result (33) gives when it is interpreted in this way. 
The potentials given by (33) are 


ky’). (33) 


1 
A(K) =4](®)| 


— f 50) exp (— 10K, exp (i(K, 


“(2 )* 
Hence, if we write ’ 
X’)= f exp (i(K, R))d*k, (35) 


the retarded and advanced potentials may be written in the form 


A(K)= (36) | 


Equating the right-hand sides of (19) and (36), we find 
5(ror) 


g(X, X’) = (37) 


Let us now find what expressions for g(X, X’) are given by (35) when Cauchy’s principal values are ‘ 
taken. { 
Consider first the integral | ' 


| 1 
kox)dky =— ( ikoxodk 


Since 
xp (ikoto)dko= exp (ike) (#20), 


exp (ikoxo)dko= +i exp (—ikxo) (x020), 


we have 


1 sin 
f (#120), (38) 


169 
tials, 
es of 
the 
ond 
) be 
(26) 
(27) 
(28) 
can 
urce 
‘the 
q 
| 1 
7 } 
| 


S. T. MA 


J = exp (i(K, X))d*k=F2m'A(X) (xoZ0), 


ky? — 


by (4). From (5), (35), and (39) we obtain 
g(X, X’)=A(R)  (ro>0), 
=0 (<0), 


o(X,X’)=0 —(ro>0), 
=-A(R) (r<0), 


for the retarded potentials and 


for the advanced potentials. 
On account of (1), formulae (37) and formulae (40), (41) agree except when r>=0. On the plane 


ro=0, (37) gives us 4(r)/r but (40) and (41) give us 0. Thus the two expressions for the function 
g(X, X’) are the same everywhere except at the point R=X—X’=0. 


4. FOURIER TRANSFORMS OF RETARDED AND ADVANCED POTENTIALS IN THE CASE OF «+0 
Let us now consider a U field satisfying the equation . 
(C+«*) U(X) (42) 


We assume, for simplicity, U and p to be scalar functions. 
On expressing the functions U(X) and p(X) as Fourier integrals in the form 


1 
J exp G(R, (43) 


1 
f exp (i(K, X))d*k, (44) 


U(X) = 


and substituting in (42), we find 
(x?-++k? — ko”) O(K) (45) 


The result of dividing out (45) by the factor (x?-+-k? — ko”) is not unique, since when an arbitrary linear 
combination of the functions 6[ko— and is multiplied by ko?) the 
product will vanish. Formulae (33) of Section 3 suggest that the Fourier transforms of the retarded 
and advanced U fields should be of the form 

By a piece of work analogous to that of Section 3 we find that (46) substituted in (43) gives, by (44), 


U(x) = f (47) 


which reduces to 


O(K) 


with 
exp (i(K,R))d*k, (48) 


(ro>0), 

d = 0 - (ro < ) ’ (#9) 
g(X, X’) =0 (ro>0), a 
=—4rD (ro <0), 


| 
170 
so that 
| 
| 


(39) 


(43) 
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46) 
4), 


47) 


48) 


19) 


50) 
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for the retarded and advanced potentials, respectively. It can be seen with the help of (10) and (11) 
that the formulae (49) and (50) agree with the results obtained by Bhabha’ for the meson field. 


5. OUTWARD AND INWARD MOVING WAVES 


The work of the last two sections was carried out on the same lines as Dirac’s treatment of scat- 
tering. Dirac obtained a general result which may be stated as follows: When |x| is large, . 


1 
exp (tk-x)f(k) | x*)§— (k?+ x?) jar 


ye (it's) (51) 


Dirac showed with the help of these formulae that by taking the combinations 


+x?) (k +x?) 
one can obtain wave functions which represent outward and inward moving particles, respectively. It 
was later shown by Mller’ in the nonrelativistic approximation that Dirac’s result holds not only for 
large distances from the source but throughout the whole space. 

It will be interesting to find out the connection between our formulae for retarded and advanced 
o'r and Dirac’s formulae for outward and inward moving waves. From (43) and (46) we have 


| 

K 


The asymptotic values of the right-hand side of (52) for large distances are, according to (51) applied 
with —x in place of x, 


exp (iL koxo= (ko? — ord, 


U(X) f (ky? — x2)! 


x 
x RoxoF (ko? — x*x)! 
+2" f exp (i[koxo — x*x) (53) 
x 
showing that the retarded and advanced potentials correspond to the — and inward moving 4 


waves, respectively, both for frequencies ko>« and for frequencies ko< 

Besides the retarded and advanced potentials, there is another type a solutions of the inhomo- 
geneous wave equations which have a simple mathematical form, namely, the potentials whose 
Fourier transforms are _ instead of by (46), by 


=4np(K) 


0(K) =41rp(K) 


x?) 


Proceeding as in deriving (53) we find as the asymptotic values of the U fields whose Fourier trans- 
Proc. Roy. Soc. A172, 384 (1939). | 


* Mgller, Zeits. f. Physik 66, 513 (1930). 
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forms are given by (54): 
{+ + (ko? — x?)! exp (i[ (Ro? — x?) (55) 


x 


Unlike the retarded and advanced potentials, this type of solutions either represents outward moving 
waves of positive frequencies together with inward moving waves of negative frequencies or represents 

‘ inward moving waves of positive frequencies together with outward moving waves of negative 
frequencies, so that waves with frequencies of opposite signs travel in opposite directions. Such U 
fields may be written in the form of (47) with 


g(X, X’) —f exp (i(K, R))d*k. (56) 
24 
These expressions reduce to 


respectively, which are Hankel functions. 
APPENDIX. PRINCIPAL VALUES OF IMPROPER INTEGRALS 


Concerning Cauchy’s principal value of an improper integral, the following theorem is frequently 
useful. Let f(x) be a function such that the integral, 


f(x) 


exists, where 


fo(x) = 
db>0. Th 
an en Se), folt) 
a—b 
For, if we write 
f(x) +f(2a—x) 


2 


As a corollary to this theorem we note that if f(x) be expressible as Taylor's series in the form 
f(x) = f(a) +(x—a) 


and 6 is a small number, then 
f(x) 
P| (60) 
a> X—a 
This corollary has been used in a recent paper by Hsiieh and the present writer.® 
I wish to express my sincere appreciation of the friendly action of the United States Government in 
sending scientific microfilms to Chinese Universities in connection with the Cultural Relations 


Program of the Department of State which has greatly facilitated this work. 
® Ma and Hsiieh, Proc. Camb. Phil. Soc. 40 (1944). 
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The Analysis of the Vibration-Rotation Band «; for C“O", and C¥O™, 
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The vibration-rotation band w; for normal and isotopic CO; has been completely remeasured 
with a 7200 line-per-inch grating and the rotation lines of both bands are completely resolved. 
Identification of the rotation lines up to J=90 in the P branch and J=56 in the R branch for 
the band in normal CO: was made, while for the isotopic band they are identified to J=94 
in the R branch and to J=44 in the P branch. The combination principle was applied to the 
new data and new values of the band centers and the rotational constants Booo, Bon, and as 
have been determined. These new values are compared with those obtained from the 


earlier data. 


I. INTRODUCTION 


N connection with the investigation of w: of 

CD;Cl at 2285 cm, it was observed that the 
rotation lines of the isotopic CO, band, with 
which it overlaps, were completely resolved. As 
it was apparent that the resolution of the lines 
was much better than in the earlier work! on 
normal and isotopic COs, it was deemed worth- 
while to reinvestigate the region of ws for both 
molecules exhaustively. 

The bands were observed on a prism-grating 
spectrometer with a two-meter-focus collimating 
mirror employing a 7200 line-per-inch replica 
grating made by R. W. Wood. Under ordinary 
conditions the amount of CO, in the atmosphere 
is so great that no transmission of radiant energy 
occurs in the most intense portions of ws for 
normal CQ,. In order to map these regions it was 
found necessary to place large trays of Ascarite in 
the spectrometer and seal all cracks in the 
enclosing box with scotch tape. This device so 
reduced the CO, content of the air within the box 
that about 40 percent of the energy was now 
transmitted in the most intense regions. 

The spectrometer slits were set to include a 
spectral interval of about 0.25 cm™, and the 
galvanometer deflections were read at frequency 
intervals of 0.07 cm~!. These deflections were 
plotted as ordinates versus the angular position 
of the grating with respect to the central image 
as abscissae. The angle scale was later converted 


into a scale of frequencies in cm~ by the usual. 


* The University of Tennessee, Knoxville, Tennessee. 

** The National University of Peking, Kunming, China. 

1P. E. Martin and E. F. Barker, Phys. Rev. 41, 291 
(1932); D. M. Cameron and H. H. Nielsen, Phys. Rev. 53, 
246 (1938); A. H. Nielsen, Phys. Rev. 53, 983 (1938). 


relation vy (cm~') = K, (cm~) csc 6, where K, is the 
instrument constant corrected to vacuum. K, 
was determined from this relation by making a 
great many observations on the angular position 
of an intense green line of the mercury spectrum 
in the sixth order. 

At least three complete sets of observations 
were made on the two bands in this region. The 
angular position of each line was determined 
from about 10 or 12 observations and in most 
cases are accurate within +0.07 cm. The ex- 
ceptions occur where a line of normal and isotopic 
CO, overlap. In the R branch of the normal CO, 
band; where the lines are closest together because 
of the convergence, lines separated by only 0.8 
cm~! were completely resolved. Where the over- 
lapping occurs, lines separated by as little as 0.2 
cm~ were consistently, though not completely, 
resolved. 


Il. EXPERIMENTAL RESULTS 
The Vibration w; for Normal CO, 


In Fig. 1 is Shown the graphical representation 
of ws for C"O"*, with galvanometer deflections 
plotted as ordinates and frequencies in cm™ as 
abscissae. It may be seen that the resolution of 
the rotational structure is complete. The faint 
lines appearing as satellites and in between the 
main absorption peaks are real, as they have been 
repeated in all the measurements. Those lines not 
attributable to C"O"*®, may be caused by ws of 
C®O!40!8, the center of which should lie in this 
region of the spectrum. Assignment of J values 
has been possible to J=90 in the P branch and 
to J=56 in the R branch. The J value given is 
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ANALYSIS OF VIBRATION-ROTATION BAND 


TABLE I. The rotation lines of ws for C%O"*,. 


Positive branch (cm~) 
No. cation v(obs.) »(calc.) 


Negative branch (cm~) 
Line Identifi- 
No. cation »(obs.) »(calc.) diff. 


2350.73 2350.72 +.01 

103 R(2) 2352.27 2352.23 + 
104 R(4) 2353.77 2353.74 + 
105 R(6) 2355.20 2355.23 — 
106 R(8) 2356.75 2356.69 +. 
107 R(10) 2358.11 2358.12 - 
108 R(12) 2359.59 2359.53 + 
109 R(14) 2360.91 2360.91 

110 R(16) 2362.27 2362.27 

111 R(18) 2363.61 2363.61 

112 R(20) 2364.94 2364.91 + 
113 R(22) 2366.24 2366.20 + 
114 R(24) 2367.49 2367.46 + 
115 R(26) 2368.75 2368.70 + 


116 R(28) 2369.90 2369.91 
117 R(30) 2371.09 2371.09 
118 R(32) 2372.24 2372.25 
119 R(34) 2373.39 2373.39 
120 R(36) 2374.51 2374.51 
121 R(38) 2375.62 2375.58 + 

122 R(40) 2376.55 2376.65 - 

123 R(42) 2377.73 2377.69 +. 
124 R(44) 2378.71 2378.70 > 

125 R(46) 2379.71 2379.68 +.03 


126 R(48) 2380.60 2380.65 — .05 

127 R(5S0) 2381.57 2381.59 — .02 

128 R(S2) 2382.47 2382.50 — .03 

129 R(54) 2383.32 2383.38 — .06 
2384.25 


99 P(2) 2348.41 2348.36 +.05 
98 P(4) 2346.78 2346.76 +.02 
97 P(6) 2345.09 2345.15 — .06 
96 Pid 2343.38 2343.50 —.12 
95 P(10) 2341.75 2341.83 — .08 
94 2340.03: 2340.13 —.10 
93 P(14 2338.37 2338.42 —.05 
92 P(16) 2336.63 2336.68 —.05 
91 P(18) 2334.90 2334.91 —.01 
90 P(20) 2333.06 2333.11 —.05 
89 P(22) 2331.26 2331.30 —.04 
88 P(24) 2329.40 2329.45 — .05 
84 P(26) 2327.49 2327.58 —.07 


l+++ +++ 


P(50) 2303.26 2303.24 + 
P(5S2) 2301.06 2301.05 


14 P(82) 2265.10 2265.21 —.11 


12 P(84) 2262.62 2262.64 — .02 
P(86) 2260.02 


P(88) 2258.39 © 
2254.68 


that for the ground state.* In the region where 
the P branch of normal CO, and the R branch of 
the isotopic molecule overlap, the positions are 
somewhat less accurately known, and identifica- 
tion is therefore a little more difficult and 
uncertain. In Table I are listed the observed 
frequencies, corrected to vacuum, of the lines in 


* The nuclear spin of the oxygen atoms is zero, and for 
this reason the part of the wave function depending on the 
nuclear spin is always symmetrical to an interchange of the 
two identical oxygen atoms. As a result of the exclusion 
or mare which requires that the total wave function shall 

ong, to a certain class, only the levels characterized by 
even J values are present in the normal state. The effect is 


~ that every other line in the spectrum is absent, i.e., only 


R(J) and P(J) with even values of J will occur. 


ws of normal CO, with the assigned J values. 
Each line in Fig. 1 is numbered and can be 
identified by referring to Table I. In adjacent 
columns of Table I are found the line frequencies 
computed with the use of constants determined 
from the data. It may be seen that the mean, 
difference is less than 0.01 cm~. The lines in 
Table I marked with an asterisk are either 
unobserved or somewhat in error because of 
overlapping with the isotopic band. 


The Vibration «; for Isotopic CO, 


Figure 2 shows the graphs of the galvanometer 
deflections versus the frequencies in cm™ for ws 


175 
78 P(30) 2323.80 2323.78 
01 
.06 ‘ 
.00 
.03 
.02 
.02 
.01 
‘ 56 P(48) 2305.41 2305.41 .00 
53 .02 
50 01 
48 .10 
46 P(56) 2296.57 2296.59 —.02 { 
44 P(58) 2294.29 2294.33 —.04 
41 P(60) 2292.02 2292.04 — .02 
38 P(62) 2289.75 2289.73 +.02 : 
36 P(64) 2287.48 2287.39 +.09 ! 
34 P(66) 2285.12 2285.03 +.09 
31 P(68) 2282.55 2282.64 —.09 
28 P(70) 2280.31 2280.22 +.09 d 
25 P(72) 2277.79 2277.79 .00 } 
23 P(74) 2275.29 2275.32 —.03 
21 P(76) 2272.84 2272.83 +.01 j 
18 P(78) 2270.56* 2270.32 +.24 
P(80) 2267.79 
] 
: 
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TABLE II. The rotation lines of ws; for C#3O'*,. 


Positive branch (cm~) 


No. v(obs.) v(calc.) diff. 


Negative branch (cm~!) 
v(obs.) v(calc.) 


2284.91 2284.92 d 
35 R(2) 2286.44 2286.44 .00 
37 R 6) 2287.84* 2287.95 —.11 
39 R(6 2289.51 2289.43 +.08 
40 R(8) 2290.81 2290.88 —.07 
42 R(10) 2292.31 2292.31 .00 
43 Ry 12) 2293.72 2293.72 00 
45 R(14) 2294.99 2295.09 —.10 
46 Riis) 2296.57 2296.45 +.12 
47 R(18 2297.88 2297.77 +.11 
48 ates 2299.02 2299.07 —.05 
49 R(22) 2300.35 2300.25 +.10 
51 rs 2301.43* 2301.60 —.17 
52 R 38 2302.35* 2302.82 —.47 
54 R(28) 2304.15* 2304.02 +.13 
55 R(30) 2305.19 2305.18 +.01 
57 R(32) 2306.40 2306.36 +.04 
58 R(34) 2307.45 2307.46 —.01 
60 R(36) 2308.54 2308.55 —.01 
61 R(38) 2309.42* 2309.63 —.11 
63 R(40) 2310.62 2310.67 —.05 
64 R(42) 2311.51* 2311.67 —.16 


2326.77 
* 


R(86 2327.64 
85 R(88) 2328.22 2328.15 +.07 
86 R(90) 2328.51 2328.46 +.06 
R(92 2328.83 
89 2329.18 


30 P(2) 2282.67 2282.56 +.11 
29 P(4) 2280.95 2280.98 — .03 
27 P(6) 2279.35 2279.38 — .03 
26 P(8 2277.72 2277.70 +.02 
24 P(10) 2275.99 2276.04 —.05 
22 P(12) 2274.23 2274.34 —.11 
20 P(14) 2272.62 2272.62 -00 
19 P(16) 2270.70* 2270.47 +.23 
17 P(18) 2269.24 2269.10 +.14 
16 P(20) 2267.30 2267.30 .00 
15 P(22) 2265.52 2265.47 +.05 
13 P(24) 2263.67 2263.62 +.05 
11 P(26) 2261.78 2261.75 +.03 
10 P(28) 2259.69 2259.85 —.16 
9 P(30) 2257.98 2257.91 +.07 
8 P(32) 2256.00 2255.95 +.05 
6 2253.93 2253.98 —.05 
5 P(36 2252.05 2251.97 +.08 
4 P(38) 2249.97 2249.95 +.02 
3 P(40) 2247.96 2247.89 +.07 
2 Pad 2245.93 2245.82 +.11 
1 2243.83 2243.70 ‘ 


of the isotopic molecule. Several lines have been 
redrawn from other data in Fig. 2 to show how 
well small details were reproduced from one set 
of observations to another. Because of the over- 
lapping of the two bands, the lines are not as 
completely resolved nor quite as accurately 
known in the R branch of this band as in the P 
branch or as in ws of the normal molecule. 
Assignment of J values to lines has, however, 
been possible to R(94) and to P(44) with a few 
exceptions where the overlapping precluded the 
observation of the individual lines. The observed 


frequencies are listed in Table I] with the ap- 
propriate J value for the lower state.? The lines 
in Fig. 2 are numbered and can be identified by 
referring-to Table II. The line frequencies com- 
puted from constants determined from the data, 
are found in adjacent columns in Table II. Be- 
cause of overlapping with the normal molecule 
band, a number of lines are either in error or 
unobserved and are marked with an asterisk. The 
average difference between the observed and 
computed frequencies is about 0.02 cm~. 


| E . 
176 
No. cation iff. 
67 R(46) 2313.52* 2313.65 13 
R(48) ° 2314.86 
69 R(50) 2315.08 2315.05 +.03 
R(52) ° 2315.96 
71 Ries} 2316.89 2316.83 +.06 
R(56 ° 2318.07 
73 R(58) 2318.90 2318.90 .00 
R(60) > 2319.70 
75 R(62) 2320.50 2320.47 +.03 
R(64) 2321.19 
R(66) ad 2321.91 
R(68) 2322.60 
77 R(70) 2322.90 2322.93 — .03 
R(72) 2323.89 
79 R(74) 2324.34* 2324.21 +.13 
80 R(76) 2324.80* 2324.99 —.19 
81 R(78) 2325.51* 2325.65 —.14 
83 R82) 3326.79 02 
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Boo, 0.3878 
Beco 0.38/0 
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Fic. 3. Plots of R(J- 1)—P(J+1) and R(J)—P(J) versus 2(2J+1) for ws of normal COs. 


Ill. DISCUSSION OF THE RESULTS 


The combination principle* has been applied to 
the bands measured, with the result that an 
accurate direct determination of the constants 
relating to the vibration ws were possible. The 
energy expressions for the rotational levels in the 
normal and first excited vibrational states are, 
respectively : 


E" /hc= 
— J" J" +1)D, 
(1) 


and 


E'/he = Eoo/he+J'(J' +1) Boos 


Booo and Boo; are constants (B=h/8x*Jc) repre- | 


senting the moments of inertia in the ground and 
first vibration states, and D is the centrifugal 
stretching constant. It may be shown for this 
molecule that, to a good approximation, D may 
be neglected. The rotation lines on the high 
frequency side of the band center form the R 


*G. Herzberg, Infrared and Raman Spectra of Polyatomic 
Molecules (Van Nostrand, New York, 1945), p. 390. 


branch and are given by: 


R(J) 
+(J+1)*(Boo1—Booo), 


while the lines on the low frequency side form the 
P branch and are given by 


P(J) J (Boor t+Booo) +J*(Boo1— Booo), 


where J is the quantum number of total angular 
momentum in the lowest state, and wo is the band 
center or the frequency of the J=0 line. By 
combining these relations properly, the following 
expressions, from which Booo and Boo: emerge, 
may be obtained: 


R(J—1) —P(J+1) =2(2J+1)Booo 
(2) 


and 
R(J) — P(J) = 2(2I+1)Boo. 


Figures 3 and 4 show graphs of these two 
equations with term differences in cm~ plotted as 
abscissae and 2(2J+1) as ordinates for both 
molecules. The slopes of these lines yield Booo and 
Boo, respectively. B,, the equilibrium value of B, 
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Fic. 4. Plots of R(J-—1)—P(J+1) and R(J)—P(J) versus 2(2J+1) for isotopic COs. 
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Fic. 5. Plot of R(J)+P(J) versus J(J+1) for normal COs. 
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Fic. 6. Plot of R(J)+P(J) versus 
J(J+1) for isotopic COs.’ 


and the equilibrium internuclear distances may 
be obtained from Booo and Boo: and the conver- 
gence factors a, a2, a3 with the relation: 


—ao(ve+ 1) —a3(vs+3). (3) 


The frequency, wo, of the J=0 line may be 
obtained from another combination relation, 
namely : 


R(J) + P(J) = 2(wo+Boo1) 
+2J(J+1)(Booi—Booo). (4) 


When R(J)+P(J) in cm is plotted versus 
J(J+1), as in Figs. 5 and 6, a straight line 
results, the intercept of which is 2(wo+Boo1) and 
the slope of which is 2(Boo:—Booo) or —2a3. The 
values of the various constants obtained from the 
present data follow in Table III. 


IV. CONCLUSION 


The Booo and Boo: found in columns 1 and 2 of 
Table III are obtained directly as the slopes of 
the curves in Figs. 3 and 4, and the differences 
between these values give a3 or Booo— Boo: in 
column 3. Another value of a3, found in column 4, 
is obtained directly from the slopes ef the curves 
given in Figs. 5 and 6. From a theoretical relation 
for as given by Nielsen‘ and utilizing values of 
the vibrational frequencies of w:, the 
potential energy constant c given by Dennison,°® 
and a value of B,, the computed value of az, 
found in column 5, was obtained. It may be seen 
that these determinations of as are in excellent 


4A. H. Nielsen, J. Chem. Phys. 11, 160 (1943). 
’D. M. Dennison, Rev. Mod. Phys. 12, 175 (1940). 
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TABLE III. Values of the constants relating to ws for C%O"*, and C#%0"*,, 


Boo (cm™) 


we (cm~!) 
band center 


0.3878 


0.3875 


C30,16 


2349.93 
2284.14 


0.0031 


agreement with each other, but do not agree with 
a3=0.0029 cm= given by Herzberg.* The combi- 
nation principle applied to Martin and Barker's! 
data on wz yields a value for Booo of 0.3903. This 
is in reasonably good agreement with the value in 
Table III. The value Booo=0.3895 cm given by 
Herzberg* appears to be low. When the combi- 
nation principle is applied to Barker and Adel’s 
data’ on v3—(v1, 2v2) it yields values for Boo: of 
Boo: = 0.3861 (10.414 band) and Boo: = 0.3862 


* G. Herzberg, reference 3, 394-395. 
7 E. F. Barker and A. Adel, Phys. Rev. 44, 185 (1933). 


cm (9.40u band), respectively. The values of 
Boo. from Barker and Adel’s data, although in 
agreement with each other, are not consistent 
with the Boo; listed in Table III, found from the 
present data on ws. 

It may also be seen from Table III that the 
constants Booo, Boo1, and a3 for normal COz are 
very nearly equal to the values of these constants 
for the isotopic molecule. This is rather good 
evidence that the equilibrium moment of inertia 
and internuclear distance and the force constants 
are the same for the two molecules. 


" 
180 
(Boor Bun) (cm) (cm™) 
rom rom caic. trom 
Booo (cm~!) Col. 1-Col. 2 5 reference (4) 
0.0032 0.0031 m7 
0.0032 0.0032 
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The Vibration-Rotation Energies of Polyatomic Molecules 
/ 
Part II. Accidental Degeneracies 


HarRALp H, NIELSEN 
Mendenhall Laboratory of Physics, Ohio State University, Columbus, Ohio 
(Received August 22, 1945) 


The first- and second-order corrections to the vibration-rotation energies of polyatomic 
molecules are dealt with in instances where two or more vibration frequencies are accidentally 
degenerate. The method of the contact transformation employed in Part I is extended and 
made applicable to the two types of first-order resonance interactions, i.e., the Fermi-Dennison 
type and the Coriolis type. The components of the energy matrix are evaluated in general, and 
examples are considered to demonstrate how the actual energies may be evaluated. 


I. INTRODUCTION 


HE vibration and rotation of a general polyatomic molecule have been dealt with quantum 

mechanically in an earlier paper, hereafter referred to as Part I,' and expressions were derived 

for the energies by a second-order perturbation calculation. The results are applicable, except in cer- 
tain anomalous cases, to specific molecular models to obtain the vibration-rotation energies. 

The perturbation method applied to this problem is related to the fact that, except for the terms 
of Coriolis interaction between degenerate vibrations and the rotation of the molecule, the terms in 
the first-order Hamiltonian, H;, cannot in general contribute to the energy before in second order of 
approximation. It is, therefore, effective to transform the Hamiltonian, H;, of the molecule by a 
contact transformation into T H T-!= Ho’ +eH,;'+éH,'--- where H,’ will contain to second order 
of approximation only the above mentioned Coriolis interaction terms. Linear combinations vf the 
zero-order wave functions may always be found so that the matrix of H,’ will have only diagonal 
elements. Using the stabilized wave functions of H,’, i.e., the linear combination of the zero-order 
wave functions which will diagonalize Ho’+«H,’, the calculation of the second-order energies is re- 
duced, in effect, to a first-order perturbation calculation. 

Denoting by 7(e) the transformation function e*S, we have to second order of approximation : 


(1) 


Transformation of leads to THT“ = H’ = + eH’ where it is readily 
shown that ; 


Ho =Ho, Hy He! — S]. (2) 


The method is very effective except in such cases where peculiar relationships exist between 
the vibration frequencies. For example, if the first-order Hamiltonian of a molecule has a term 
Res'sQeoe'e’Qe’’e’’ in it, the transformed second Hamiltonian, H,’, will contain a contribution inversely 
proportional to 


+ (A. dart) (Ast (A. et = J, 


where A, =47°c’w,?, and the anharmonic coefficients x,,, Xs" etc., in the expression for the vibration 
energy of the molecule will contain terms with the same denominator (see for example, Eq. (30), 
Part I). If accidentally should be very nearly equal to (i.€., in which case 
resonance is said to occur, the quantities x,, and x,., may become indefinitely large and the method 
_ described above will fail. A similar instance may arise when the first-order Hamiltonian, H;, of a 
polyatomic molecule contains a Coriolis interaction term 


1H. H. Nielsen, Phys. Rev. 60, 794 (1941). 
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The transformed second-order Hamiltonian, H;2’, will contain terms corresponding to those with the 
denominator A,—),’. These in turn give rise to correction terms to the effective reciprocals of inertia 
with the same denominator. As before, if \, should accidentally be very nearly equal to d,- so that 
resonance sets in, the correction terms may become indefinitely large, and the method described before 
becomes inadequate. We shall hereafter consider how the above method may be modified to take 
account satisfactorily of the contribution to the energies by such terms when resonance between 


frequencies takes place. 


Il. ANHARMONIC RESONANCE INTERACTIONS 


We shall first consider the case where the resonance interaction results because of the presence of 
terms in the anharmonic portion of the potential energy which are cubic in the coordinates. The 
general anharmonic potential energy term which is cubic in the coordinates may be written* 


atts (3) 


This term may be removed from the first-order Hamiltonian by transforming H with the transforma- 
tion function® 


The transformation produces, however, a term in the second-order Hamiltonian with the denominator 
(As*-++Ay-#—A,--#). When, for example, the combination frequency A,!+),-! (or overtone frequency when 
s=s’) and the fundamental \,--! are very nearly alike, resonance sets in and the term Ackss’s'QseQs'e’Qs''e"” 
may evidently contribute to the energy of the molecule an amount which is of first-order magnitude. 
It therefore appears that it is not legitimate to remove the term [cRss14""QseQs'e’s'’e’ from the first- 
order Hamiltonian. A closer inspection reveals, however, that it is only the Use", | Use 1, 
Vee? 1, Ugg" 1) matrix elements of this potential energy which offers difficulties. It therefore be- 
comes apparent that if a transformation function S can be found such that after the transformation of 
Hi has taken place will contain only the elements (dso, Use", | Vso 1, 1, F 1), these 
having, of course, the same values as the corresponding elements of the term hcRgs's'QseQs’e’Qs’te"’, then 
a transformation method may still be employed with advantage. It will be seen that the term, 


has the same (0.6, | 1, 1, Vee? 1) matrix elements as (3), but that all its other 
matrix elements are equal to zero. Since the function (4) will remove the term (3) entirely from Hy’, 
we may expect to obtain a transformation function appropriate for us by adding (4) to a function s 
which will put the expression (5) back into the transformed iy’. A little experimentation shows that 


the function 
X { — Perret /It) + + /h)} 
Aad — (Nat — Aart) (At Aart) (6) 


will accomplish this. The transformation function S* suitable for our problem will be S*=S+s and 


2 The coordinates 9,, used here are dimensionless so that if ky." is expressed in cm~, the term hck,.’,’« will have the 
dimensions of energy. The constants k,,.’.” used here will be found to be equivalent to the k..’.’ used in Part I, multiplied 


by the constant 
3 W. H. Shaffer, H. H. Nielsen, and L. H. Thomas, Phys. Rev. 56, 895 (1939). S. Silver and W. H. Shaffer, J. Chem. 


Phys. 9, 599 (194i). 
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is equal to 
St= / 4) { (As — t+ t+ h‘) — 
X } { (Ast (Ast — Aart) (Ast (7) 


It should be noted that S* does not contain the denominator (A,4+A,/?—A,4). 
We have, then, for the zero-order Hamiltonian 


Ho= (1/2) + (Py?/Lyy) + +(h/2) Ds De (8) 

and for the transformed first-order Hamiltonian, H,’, the following: 
H = 4) { + Parte’ / (PeoPs'e'/ — } ° (9) 
The transformed second-order Hamiltonian will also differ from the relation given in Eq. (21c), 


Part I,‘ but only insofar as the terms containing a as a coefficient are concerned. We restate, 
therefore, only this term. It is: 


(i/2) (S*H, — H,S*) = 2 + 2A t+ 


x { 2 (Ast + + Agr?) (As } (10) 


To arrive at the energies we proceed now by assuming that A,!+A,-?~),--? so that we may write 
Ay-t=A,-?+A,!+6 where 6 is small. When we substitute this into the portion of H» which concerns 
these frequencies we obtain ' 


H'= (6h /2) gure’) { + Porras /h? 
— — (12) 
If now the frequencies are non-degenerate so that we may set ¢=o’=1, we obtain at once that 
(Eo/he) = ws (ve +1) (ve +00+1). (13) 


The first term in H;’ has only elements diagonal in the quantum numbers and the second, as we have 
seen, has only the matrix elements (v,, +1, 1). These elements are the following 


Ver, Vere | /he| ve, Ver, Vere) = (Ver A=(5/4rc), (14a) 
Ver, Vere | /he|vet+1, ve +1, —1) = (ve +1, ve +1, —1| /he| v6, ve, 
= ((Ve+1)/2) +1)/2) (00/2). (14b) 


The matrix components of the transformed second-order Hamiltonian which are of importance 
are those diagonal in the quantum numbers. In our case they are the same as those given in Part I 


‘To arrive at the modified expression for H,’, it is necessary to recompute the expression (i/2){S*(H:’+H,) 
It is quickly seen that S* has no matrix elements of the type +1, while Hy’ 
was so determined that it contains only elements of this type. In evaluating the energies to second-order of approximation, 
we shall be interested only in the elements of H2’ diagonal in the quantum numbers »,, 2’, 2%, but it is evident that 
(S*H,'— H,'S*) will have none such. It is necessary, therefore, only to consider (S*H;—H,S*). 
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except insofar as the terms having the coefficient kuv'ae are concerned. These occur in only xs. In 


our case the term multiplied by Reve will be 


1 1 1 


+ 
Ws + We’ Ws — Ws’ — Ws’ 


The expression (15) will be seen to be identical to the relation occurring in x,. as defined in Part | 
except that the term having the denominator (w,+wj—w,') is absent.® 

When the frequency A,! is equal to the frequency d,-4, we have the interesting case which is generally 
referred to as Fermi-Dennison* resonance. When the frequencies \,' and ),--? are non-degenerate so 
that o=o’’=1, the term in the first-order Hamiltonian which allows the resonance interaction to 
occur, may be written hcksss’’g."ge. The appropriate transformation of this term may be achieved by 
using the function (7) where s=s’, o=o’, and where the order in which the operators occur is pre- 
served. In this case we set A,’ =2A,!+-4, 6 being small, and we have for the zero-order Hamiltonian 


which relates to and 
Hy= (Wid. /2){ (16) 


and for the transformed first-order Hamiltonian 
Hy! = (5h /2) (parr? /h? Gurr?) + + Parr /h? — (p.?/h? — qu?) (17) 


For the zero-order energies we obtain at once 
(Eo/he) =w,(v.+ 20. +3). (18) 


As before the first term in H,’ has only elements diagonal in the quantum numbers, while the second 
_ has the elements 1) only. These are 


(vs, ver | Hi’ /he| Vs, Ver) =A(ve +4), (A=8/4rc), (19a) 
(vs, ve | Hi’ /he|v.+2, ve —1) = (0.42, | Hy’ /he| vs, 
= (04/2)  (19b) 


The diagonal elements of H,’ are again the only ones which are of interest. Their contributions 
to the energy can be shown to be identically the same as those give in Part I except insofar as the 


contributions which have the coefficient k,...- are concerned. Such terms occur only in X,, and X44". 


As in the earlier example one may show that the contributions which have-the coefficient Rese’ 
remain the same as those given in Part I except that the term which has the denominator (2, —w,’) 
will be absent. 

When the frequency A,! is twofold degenerate, as may be the case in linear and axially symmetric 
polyatomic molecules, the index ¢ may take the values ¢=1 and o=2. In such cases resonance may 
occur between the overtone 2w, and the fundamental frequency w,- if the first-order potential energy 


contains the following component ; 
oss’ (Qo, 27) Qs". (20) 


Here it is only the elements (0,13 Us, 23 | Ve, F 1) AN (Os, 15 | Vs, 15 Vs, 22, Ve 
which cause trouble when resonance between 2w, and w, occurs. The terms in (20) are properly 
transformed by use of (7) exactly as in the preceding example. As before we proceed by setting 


5 The component of x4." which has as a coefficient is stated in Part I as follows: / tees) 


XK (wet — (oe — (Wop — We’ — This may be factored into (1/4) 


1 


* E. Fermi, Zeits. f. Physik 71, 250 (1931). D. M. Dennison, Phys. Rev. 41, 304 (1932). 
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hy? =2A,4+6 and obtain for the zero-order Hamiltonian the following : 
Hy = { (ps, + Pa, 1° +90, 2") +2 /h? }, (21) 
and the transformed first-order Hamiltonian will be: | 
Hy! = (54/2) (Pare? /h? + /4) 19, 149s, 1+ Pe, 29,2 
21,2) Par /h? — (Da, + Pa, 2?/h? — 1? — Ga, 2") (22) 


It is convenient to replace the coordinates g,1 and g,,2 and the conjugate momenta p, and 
ps2 by their equivalents in cylindrical coordinates. These are COS Xe, SIN Xe," 
= (—ih)[cos x,0/dr,— (sin and p,2=(—ih)[sin x.0/0r.+(cos The eigen- 
values of (21) are known to be: 


(Eo/he) = (23) 
The matrix components of H,’ rewritten in the same coordinates are known to be the following :” 
(vs, | Hy’ /he | Vs, =A(v,+4), (A= 5/4xc) (24a) 


Vere | /he| v.42, Le, ver —1) = +2, Le, Ver —1| Hy’ /he|v., Liver) 


(—Reser/2)[ (02 +2)?—1,? 

where /, takes the values v,, v,—2, ---1, or 0. 
As in the preceding example the elements of the transformed second-order Hamiltonian which 
are diagonal in the quantum numbers are the same as those in Part I, except insofar as those terms 


which have as coefficients are concerned. These elements contribute only to ANG 
and the terms which have as a coefficient are altered from (— 40411) (802 — / (402 — 
to — ] in Xsey from — (eas? — to (— in 


tur, and from to in These will be seen to be 
the same as those occurring in the expressions for 4s, Xs.’", and Xi,/, given in Part I, except that the 
part containing (2w,—w,’,) as a denominator is absent. 

The actual energy values may be obtained by diagonalizing the matrix H, and this is readily ac- 
complished by setting its secular determinant (i.e., |(7|H|r’)—é,,-| where + comprises all the 
quantum numbers, vibrational and rotational) equal to zero and solving for the roots. For simplicity 
we confine ourselves to the energies of the non-rotating molecule. 

Study the first example considered. The elements of H are diagonal in all the vibration quantum 
numbers except v,, v,", and v,. The zero-order energy (13) no longer depends upon the quantum 
numbers 2,, v,, and v,- independently, but rather upon v,+2,, and v,-+9,--. There exists, therefore, a 
degeneracy and the degree of degeneracy is the number of ways v,+0, and v.+2,-- can be made to 
add up to a given value. For example, consider that v,+v, and v,-+2,- are made equal to one. This 
can happen if v, =v, =1 and =0 or if v, =v, =0 and v,,=1. In both cases (Eo/hc) = and 
the level is twofold degenerate, and the matrix elements relating to this level form a step matrix of 
two rows and columns grouped about the principal diagonal. This is illustrated in (25): 


t.=v7=1, v=1 
v,=0, =1 
0 2(wet+ws’) +A/2 Resa 
(25) 
v.=v,=0 
/(8)3 2 (wet we) +34/2 


? These matrix components have been computed by Professor W. H. Shaffer of this University who mg made his 
computations available to the author. The author wishes to express his appreciation to Professor Shafter 


for his kindness. 
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Similarly, in the second example the elements of H are diagonal in the vibration quantum numbers 
except v, and »,. There will, as before, be a degeneracy since the zero-order energies (18) do not 
depend upon the quantum numbers, », and »,-, independently, but upon v,+2v,. The degree of the 
degeneracy is determined by the number of ways in which v,+ 22, can be made to add up to a given 
value. Suppose w, is a linear harmonic frequency and that we let 2v,--++-v, be equal to four. This may 
take place in the following ways: =0, v, = 4; =1, v, =2; and =2, v,=0. The level 
is, therefore, threefold degenerate, and the matrix elements relating to this level form a step matrix of 
three rows and columns grouped about the principal diagonal. This step matrix is illustrated in (26) ; 


v,=0, v,=4 ve =1, v,=2 =2, v,=0 


=0 

| | Chore /2)(6)3 0 
Vy =2 


The example where w, is a two-dimensionally isotropic frequency has been examined by Dennison* 
and will not be discussed here. 

The matrix of the vibration energy will, therefore, to this approximation contain elements along the 
principal diagonal only, except insofar as 2,, v,-, and v, are concerned. The diagonalization may then 
be accomplished by diagonalizing each substep independently. When the roots have been evaluated, 
the stabilized wave functions for the component states may readily be obtained since these will be 
linear combinations of the wave functions of the unperturbed states, the coefficients multiplying each 
term being the normalized cofactors of the secular determinant. 


Ill. THE CORIOLIS TYPE RESONANCE INTERACTION 


We shall investigate, in this section, a perturbation which arises when two mutually perpendicular 
vibrations interact with the rotational motion of the molecule through the Coriolis operator : 


IL zz) + (nese P »/ I vw) + (¢ :/ I es) J. (27) 


The quantities £ses'6’, Nsos’e’y AN {seee Which occur in (27) are the Coriolis coupling coefficients and 
depend in an involved manner upon the nature of the normal coordinates associated with the fre- 
quencies w, and w,. The term (27) may be removed from the first-order Hamiltonian by transforming 


H using the transformation function: 
S= — + 2st ret As’) } 
x { I ss) + (Meese P »/ vy) + (¢ I ss) } (28) 


The transformation of H by S gives rise to terms in the second-order Hamiltonian with the denomi- 
nator (As—Ag4’), (i.e., ws — ws’). When X,! is very nearly equal to \,-!, resonance sets in and we have much 
the same difficulty as that which arose in the preceding section. In this example it is the 
(Vee, Vere 1, 1) matrix elements of the operator (27) which present difficulty. We may 
proceed as in Section II to transform H by a contact transformation such that it will have only the 
elements (vse, Us'e’ | Vee 1, Vere’ 1), these being the same identically as the corresponding elements of 


*D. M. Dennison, Rev. Mod. Phys. 12, 175 (1940). 
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the operator (27). It will be noted that 


+ y/ I vy) + I os], (29) 


satisfies the above requirements. Since S will remove the operator (27) completely from Hy’ we may 
expect to derive the function S* required by adding to Sa function s which will put the expression (29) 
back into the transformed first-order Hamiltonian H;’. It may be verified that the correct form for s 


will be 
Py /Tyy) + (30) 
Evidently the transformation function S* =S+s suitable for our problem will be: | 
+ P y/Tyy) + (31) 


It does not contain the denominator (A,!—),/#). 

Consider now the two frequencies w, and w,, which are assumed to arise from two linear harmonic 
oscillations which are perpendicular to each other. We proceed as in Section II by assuming that 
h,'=A,/? so that we may write \,-?=A,!+4, 6 being a small quantity. The zero-order Hamiltonian 
pertaining to these two frequencies may be written: 


Hy = (hd/2)[ (ps? /h? + (per?/h? J, (32) 
(Eo/he) = (33) 
The first-order transformed Hamiltonian will be: 
Hy! = (5h /2) (pe? /h?+qe*) (34) 
Hi’ has the following non-vanishing matrix components: 
(+++) Vere K| Hy! /he| ++ K| Hy! /he| + +09, K) 
(35b) 


Hy has the eigenvalue 


Vert Ay! /he| ve —1---K+1) 
me —K(K41) PL (Ad 1) (35c) 


The transformed second-order Hamiltonian H,’ will, of course, also be altered, but only insofar as 
the coefficients of the operators P,’, P,?, and P? are concerned which are multiplied by &..:*, n..*, and 
fs. When the contact transformation of H has been completed, the above coefficients will be replaced 


by 
As) 1+3(A./ ds) h*) } { z/ I ss) + »/ I + (¢ wP 2/ I ss) } (36) 


To this approximation H,’ will contribute only to the diagonal elements of the matrix H, and these 
- will be exactly the second-order corrections to the reciprocals of inertia given in Part I except that 
the terms containing \,/(A.— As’) which occur in the rotational constants ANd Zee Will be absent 
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and in their place will be found, respectively : 
neat? yy (Ast — Aart) J, 


/2T (Aadar)*)[ (Ast — / (Ast + 
(37) | 
The energy values themselves are obtained as 
before by solving for the roots of the secular deter- : 
minant of the matrix which has the elements (32), 
(35), and (37). To second order of approximation 
all the terms are diagonal in the vibration quan- 
tum numbers », and »,. The zero-order energies 
are degenerate in the frequencies w, and wy, since 
— they depend not upon », and »v, independently, 
but upon v,+0,. The degree of the degeneracy 
depends upon the number of ways v,+v, can be 
made to add up to the same number. The ele- 
ments associated with the resonating components 
form a sub-matrix of the matrix of H. This sub- 
matrix will contain as many rows and columns 
as the degree of the degeneracy, its elements 
ole being grouped about the principal diagonal of 
the matrix in a little box. Diagonalization of the 
matrix is accomplished simply by diagonalizing 
these sub-matrices by themselves. Once the roots 
of the secular determinant are known, the stabi- — 
lized wave functions are immediately available 
since they are linear combinations of the zero- 
order wave functions of the component levels, 
the coefficients of these being the normalized first 
minors of the secular determinant. 
We shall examine the case where the frequen- 
cies w, and w, lie in the xy plane of the molecule 
} as is the case in the planar ZX Y2 molecular 
model studied by Silver. This is tantamount 
to setting &* =. =0 in which case the matrix 
elements (v,, 04°, K |v,+1, 1, K+1) will vanish. 
If we confine ourselves to the fundamental fre- 
quencies w, and w» (i.e., =1, =0, and v,=0, 
= 1), we have specifically the case described by 
Silver.* We shall set up the secular determinant 
of the sub-matrix associated with these frequen- 
cies for the value of J=1, J being the quantum 
number of total angular momentum. We follow 
Silver and use as the basic rotational wave func- f 
tions (})'[y(K =1)+¥(K=—1)] and ¥(K=0) 
instead of the functions ¥(K = 1), ¥(K = —1) and 
¥(K=0). We have then Eq. (38) as given:in rs 


*S. Silver, J. Chem. Phys. 10, 565 (1942). 81 
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adjacent column where the R; 
are the values of R; defined in 
Eq. (27), Part I, where v,=1 and 
y, =0 and the R;,’ are the values 
when v, =0 and v,, =1. The roots 
of (38) will be seen to be those 
given by Silver® for the state 
J=1. 

The selection rules for the 
quantum number K are altered 
by this type of resonance inter- 
action. This question has been 
dealt with, elsewhere!® and need 
not be repeated here. 


4(Root+ Root 
—Ra—R»— 

4(Roa+ Reo 
+2¢Z.—e 


4(Roat+ Rao 
—Ra—Ra—Re) 

4(Roat+ Rae 
+2¢Z.—e 


Another example, which may 2 
be regarded as a generalization of + 
the previous example, but which | 
is sufficiently different from the 
outs foregoing to warrant separate 
umes mention, is the one where there 
nani are three resonating frequencies, 
ad We, and wy (i.e, 
f the and 5 and being 
izing small quantities), arising from 
seat three linear harmonic oscillations 
which are mutually perpendicu- 
lable lar to each other. We shall con- 5 5 

fine ourselves to the situation me me 
zero- + | 
vels, where the vibrations associated 
with the frequencies w,, w,., and = 

w, are respectively along the ~ 
men body fixed axes z, x, and y. This = alt 
is equivalent to declaring that 
only the Coriolis coupling factors, 
anil and are different + + 
from zero. In this case there will e+ | 
nich: be matrix elements of all the types = 
Y™ given in (35) which are non-van- 
=0 ishing. We shall here consider only Hi 
4 by the fundamental frequencies, w., 
and (i.€., 0, =1, =0; 
uen- = 0, = 1; =v =0, 
wane v=1), and set up the secular 


determinant for the sub-matrix 
sal associated with these frequencies 
=0) for the value of J=1, J being the 
| ol quantum number of total angular 
momentum. We follow a method 


1H. H. Nielsen, J. Chem. Phys. 5, 
818 (1937). 
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suggested by the work of Shaffer, Nielsen, and Thomas’ and use as the basic wave functions the 
1 0 
— F R(K), G R(K), and H R(K), F and H being, (4)! 1) and G 


where =0, =1, =0) and ve =0, =0). In the 


above atti R(K) is the vila wave function which for a given J value is characterized by 
the quantum number K. We have then the rather complicated secular determinant" given in Eq. 
(39) on page 189, where 


B= (1/4) {LV =1) + =1)]-[X 0" =1) =1)]}, 


and the R,’s are defined as before. 

When the molecule is a symmetric molecule, the reciprocals of inertia X, and Y, are alike, the 
frequencies w, and w,” are the same (i.e., As4+=A,--#), and &, and 74", become equal to each other. The 
determinant (39) is then much simplified since 8 and y reduce to zero and the R;, become respectively 
equal to the Ri». It will then be seen that the secular determinant factors into five sub-determinants, 
two identical ones with a single element, two other identical steps with two rows and columns each, 
and another with three rows and columns. The general secular determinant for a given value of J 
will have 3(2/+1) rows and columns, but will again factor into sub-determinants. There will always 
be two identical steps with a single element each, two identical steps each with two rows and columns, 
and the rest will have three rows and columns each. The general form for the sub-determinant will be: 


(40) 


F R(K+1) G R(K) H R(K—-1) 
Ro+(K+1)*Rz . 
FR(K+1)| +(K+1)'Rs 0 
G R(K) Roet+K?R2- =0. (41) 


f-K(K-1) Ro+(K—1)?R2 
H R(K—-1) 0 -| +(K—1)*R; 
+2¢(K—1)Z,.—e 


There will be an odd number of steps with three rows and columns each, and these will be identical in 
pairs except for the step where K =0. In the limiting case of molecules like CH, where X,= Y,=Z,, 
the frequencies w,, w,, and w, form a triply degenerate frequency. It may be shown that £&y'5= 4 
={,,=¢. If the centrifugal distortion term R3 is neglected in the sub-determinant (41), it yields 
(e—Ro)/he=2S¢B., (e—Ro)/he= —2¢B., and (e—Ro)/hce =-—2(J+1)¢B, as the values for the roots. 
These are the well-known values given by several authors.” 

It is not practicable here to discuss the selection rules which result because of the resonance since 


they become quite complicated. One may visualize them somewhat in the following manner. The 


The elements (F R(1) (H R(1)), (F R(—1)|HR(1)), and (H R(1)| F R(—1)) are not exactly 
zero as indicated in (39). are of the o ~ Ni differences in the anharmonic corrections to the reciprocals of inertia 
and will, therefore, not thane ute to the tomy | of the molecule in this approximation. 

2M. Johnston and D. M. Dennison, Ph 


ev. 48, 868 (1935). H. A. Jahn, Proc. Roy. Soc. A168, 469 (1938). 
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frequencies w, and w, were assumed to lie in the xy plane of the molecule while the frequency w, was 
taken to be along the body fixed axis z. When the molecule is a symmetric molecule and w, = w,", we 
have a twofold degenerate frequency in a plane perpendicular to the axis of symmetry. Such a fre- 
quency is termed a perpendicular vibration and the selection rules in the unperturbed case are 
AJ=+1, 0; AK=+1. The frequency w, is non-degenerate and along the axis of symmetry. It is 
known as a parallel vibration and in the unperturbed case, the selection rules are AJ = +1, 0; AK =0. 
When resonance sets in there will be a mixing of the wave functions, i.e., the wave functions of the 
component levels become linear combinations of the unperturbed wave functions. Thus, in the case 
discussed here, the actual wave functions will be linear combinations of F R(K+1), G R(K), and 
H R(K-—1). This will mean that when a transition takes place from the normal state to the component 
level which corresponds more nearly to the unperturbed state, w,, the resulting spectrum will, to be 
sure, have the characteristics of a parallel type band, but at the same time it will also possess charac- 
teristics of a perpendicular band. In the same manner when a transition occurs from the normal state 
to the component level w,(=w,’-) the spectrum will have characteristics not only of a perpendicular 
band, but also to some extent the characteristics of a parallel type band. 

By way of conclusion it is perhaps worth while to point out the similarities and the differences 
between the method here employed to obtain the vibration-rotation energies when vibration fre- 
quencies are accidentally degenerate and the method used by Dennison" in his study of the carbon 
dioxide molecule. The method of Dennison is to evaluate the energy matrix by using the zero-order 
wave functions, i.e., the solutions to the zero-order wave equation. This matrix will, to the desired 
approximation, have diagonal elements containing zero- and second-order contributions and non- 
diagonal elements of first- and second-order magnitude. This matrix is then transformed into another 
matrix which has no off diagonal elements of less than second-order magnitude. To second order of 
approximation the elements along the principal diagonal may be regarded as the energies of the 
molecule. When two vibration frequencies, w, and rw, are close together, the method may fail since 
the transformed matrix will have second-order correction terms on the principal diagonal which have 
the denominators w,—7w,’. In such instances where w, ~ Tw, the elements which cause trouble are the 
same as those referred to in earlier sections. These are computed on the basis that the energies E, 
and E, are degenerate, and the corresponding elements will be contained in sub-matrices grouped 
about the principal diagonal. The remaining elements are computed as before, the elements along the 
principal diagonal being the same, according to Dennison, except that all terms containing the 
vanishing denominators will be absent. 

In our method it is the Hamiltonian which is transformed in such a manner that the first-order 
Hamiltonian, //;’, vanishes to second order of approximation so that the second-order corrections 
to the energy are simply the diagonal elements of H,’, and consequently the energies of the molecule 
may be taken to be just the diagonal elements of H)+eH,’. When resonance occurs the transforma- 
tion is of such a character that the first-order Hamiltonian, H;’, while it does not vanish, contains 
only terms which would give rise to energy corrections that have vanishing denominators. As in 
Dennison’s method their contributions are determined by degenerate perturbation theory, and the 
elements may be arranged about the principal diagonal in little boxes. It is readily verified by this 
method that the diagonal elements are exactly the same as when no resonance exists except that they 
contain no contribution with vanishing denominators. 

The author wishes to acknowledge his indebtedness to Dr. Samuel Silver of the Radiation Labora- 
tory at the Massachusetts Institute of Technology who has read this manuscript and verified many 
of the relations included in it. 


4D. M. Dennison, Rev. Mod. Phys. 12, 175 (1940). 
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